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Abstract 

A numerical method is developed for solving periodic, three-dimensional, vor- 
tical flows around lifting airfoils in subsonic flow. The first-order method that is 
presented fully accounts for the distortion effects of the nonuniform mean flow on 
the convected upstream vortical disturbances. The unsteady velocity is split into a 
vortical component which is a known function of the upstream flow conditions and 
the Lagrangian coordinates of the mean flow, and an irrotational field whose poten- 
tial satisfies a nonconstant-coefficient, inhomogeneous, convective wave equation. 
Using an elliptic coordinate transformation, the unsteady boundary value problem 
is solved in the- frequency domain on grids which are determined as a function of 

the Mach number and reduced frequency. 

The numerical scheme is validated through extensive comparisons with known 
solutions to unsteady vortical flow problems. In general, it is seen that the agree- 
ment between the numerical and analytical results is very good for reduced frequen- 
cies ranging from 0 to 4, and for Mach numbers ranging from .1 to .8. 

Numerical results are also presented for a wide variety of flow configurations for 
the purpose of determining the effects of airfoil thickness, angle of attack, camber, 
and Mach number on the unsteady lift and moment of airfoils subjected to periodic 
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vortical gusts. It is seen that each of these parameters can have a significant effect 
on the unsteady airfoil response to the incident disturbances, and that the effect de- 
pends strongly upon the reduced frequency and the dimensionality of the gust. For a 
one-dimensional (transverse) or two-dimensional (transverse and longitudinal) gust, 
the results indicate that airfoil thickness increases the unsteady lift and moment at 
the low reduced frequencies but decreases it at the high reduced frequencies. It is 
seen that mean airfoil loading leads to a significant reduction in the unsteady lift 
for the low reduced frequencies for both the 2-D and 3-D gust cases, but has little 
effect for the 1-D case. Finally, the results show that an increase in airfoil Mach 
number leads to a significant increase in the unsteady lift and moment for the low 
reduced frequencies, but a significant decrease for the high reduced frequencies. 
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CHAPTER I 


INTRODUCTION 


Most Sows encountered in aerodynamics are high speed flows where the 
Reynolds number is large and the effects of viscosity are confined to small regions 
such as boundary layers and wakes. Because major portions of these flow fields are 
essentially inviscid and irrotational, potential flow theory has been used extensively 
by aerodynamicists in the analysis of flows about streamlined bodies. Today steady 
potential flow solvers are widely used in the design of aircraft wings, turbomachinery 
blades, and helicopter rotors. 

In many real flow applications, however, the flow is not steady but unsteady. 
Frequently the unsteadiness in the flow is due to the occurrence of upstream vortical 
disturbances that are convected downstream and induce an unsteady flow field as 
they interact with the body. For an aircraft wing, such upstream flow distortion can 
be caused by atmospheric turbulence. For propeller and turbomachinery blades, the 
vortical disturbances may be caused by the viscous wakes of an upstream rotor or 

stator, installation effects, or upstream turbulence. 

When viewed from the blade frame of reference, the upstream vortical distur- 
bances will appear as propagating vorticity waves that axe called gusts. There are 
a number of undesirable effects that can be associated with such vortical gusts. 
They will, for example, induce unsteady forces on the airfoil surface which can 
cause forced vibrations and radiate noise into the far field. In some instances, the 
impinging gusts may cause flow separation and loss of aerodynamic performance. 
For rotating blades, the fundamental frequency of the upstream disturbances will 
equal the blade passing frequency. If the frequency of the aerodynamic excitation 
equals a natural frequency of the rotating blades and the amplitude is sufficient, 

then catastrophic structural failure may result. 

Another possible source of unsteadiness in the flow is the unsteady motion of 
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the airfoils or blades themselves. Such unsteady structural motion can be caused by 
structure-borne vibrations as well as the flow-induced oscillations described above. 
There can also be unsteady interactions between the airfoil motion and the incident 

disturbances which can dampen or increase the magnitude of the airfoil unsteady 
motion. 

Because of the undesirable effects associated with these unsteady flows, there is 
considerable interest in controlling and understanding the aerodynamic excitations 
which can cause such unsteady blade motion. 

The earliest attempts to analyze these kinds of unsteady flow fields goes back 
to the 1920’s and 30’s when aeronautical engineers sought to deal with vibration 
problems that arose due to higher flight speeds. The early work dealt with flat 
plate airfoils at zero mean incidence in incompressible flow. By considering small 
disturbances to a uniform mean flow, it was possible to linearize about the mean 
flow state and uncouple the unsteady part of the flow from the underlying steady 
flow. The basic mathematical tools used by researchers at that time were primar- 
ily conformal mapping techniques and the simple concepts of circulation theory. 
Theodorsen 1 was the first to present a complete analysis for the problem of a flat 
plate airfoil undergoing sinusoidal, oscillatory motion in an incompressible flow. 

The analysis of vortical flows about thin airfoils in incompressible flow was 
pursued along similar lines. The standard mathematical treatment consisted of de- 
composing the unsteady velocity into solenoidal and irrotational components. By 
splitting the unsteady velocity in this way, the mathematical formulation of the 
gust response problem becomes similar to that of the oscillating airfoil problem. 
The solenoidal part of the unsteady velocity represents a purely convected vorticity 
wave whose mathematical expression is a known function of the upstream vortical 
disturbances. There is no pressure associated with the vortical velocity. The irro- 
tational component of the unsteady velocity is expressed as the gradient of a scalar 
potential, and is coupled to the vortical component only through the airfoil bound- 
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ary condition. The pressure is determined entirely from the unsteady potential, 
which for incompressible flows is governed by Laplace’s equation. Sears 2 derived an 
analytical solution for the unsteady lift on a rigid flat plate airfoil moving through 
a sinusoidal vortical gust pattern in an incompressible flow. 

When the mean flow is compressible, the unsteady potential satisfies a homo- 
geneous, constant coefficient, convective wave equation. In this case the mathemat- 
ical problem can be formulated in terms of an integral equation known as Possio’s 
equation . 3 Numerical results axe usually obtained using collocation techniques. So- 
lutions of periodic vortical flows around flat plate airfoils in compressible flow can 
be found in the recent review paper of Atassi . 4 

For flows around real airfoils with thickness, camber, or angle of attack, the 
mean flow is no longer uniform. Later work in the analysis of unsteady, vortical 
flows about airfoils dealt with the problem of mean flow nonuniformities and its 
effect on the airfoil unsteady response. Horlock 5 used a heuristic approach to par- 
tially account for the effects of airfoil incidence on the unsteady lift of airfoils in 
incompressible flow. Naumann and Yeh 6 used a similar approach to account for the 
effects of airfoil camber. However, both of these treatments neglected the coupling 
between the unsteady velocity and the mean potential flow about the airfoil, and 
only took into account the modified boundary condition at the airfoil surface. 

Goldstein and Atassi 7 were the first to present a systematic theory that fully ac- 
counts for the effects of the nonuniform mean flow on the airfoil unsteady response. 
Their analysis of two-dimensional, periodic vortical flows around thin airfoils with 
small camber and angle of attack showed that the vortical waves are nonlinearly 
distorted as they are convected by the mean potential flow about the airfoil. Their 
results showed that the mean flow distortion has a strong effect on the airfoil un- 
steady response, and that the wave number in the normal direction has a significant 
effect on the unsteady lift for these kinds of distorted flows. 

In a subsequent paper 8 , Atassi showed that for thin airfoils in an incompress- 
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ible flow, it is possible to obtain the airfoil unsteady lift due to two-dimensional, 
periodic vortical gusts by linear superposition to the Sears lift of three independent 
components accounting separately for the effects of airfoil thickness, camber, and 
angle of attack. Explicit analytical formulas are presented in References 7 and 8 
for the unsteady lift due to two-dimensional, sinusoidal, vortical gusts about thin 
airfoils with small mean loading. 

By fully accounting for the effects of the nonuniform mean flow, Goldstein and 
Atassi 7 have introduced the concept of the “rapid distortion theory” of turbulence 
in unsteady aerodynamics. The rapid distortion approximation was first used by 
Prandtl 9 , Taylor 10 , Ribner and Tucker 11 , and Batchelor and Proudman 12 to study 
the change in turbulent structure in a contracting stream. It is based on the ap- 
proximation that in high speed nonuniform flows, the distortion of large eddies is 
mainly caused by mean flow variation. This leads to neglecting inter-eddy inter- 
action and subsequently to the linearization of the unsteady gover ning equations 
about the nonuniform mean flow. 

While the theoretical treatments described above have afforded significant in- 
sight in unsteady aerodynamics, they are limited to incompressible flows about 
airfoils with small mean loading and thickness, or to compressible flows around un- 
loaded, flat plate airfoils with zero thickness. In addition, the vortical disturbances 
are in general limited to the two-dimensional case, so that spanwise variation of 
the unsteady flow is not included, except for flat plate results such as are found in 
Atassi 4 and in Graham. 13 

For most flows encountered in applications, however, one deals with heavily 
loaded airfoils at high Mach number and with three-dimensional upstream gusts. 
This is particularly true for turbomachinery and propeller blades, where recent 
technological trends have led to the development of highly swept, heavily loaded 
blades that operate in the transonic flow regime. Because of the high sweep, three- 
dimensional effects are important. In addition, compressibility effects are significant 
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due to the high relative Mach numbers of the rotating blades. Due to the complica- 
tions of compressibility, heavy loading, and three-dimensional effects, the analysis 
of these kinds of unsteady flows is clearly beyond the realm of analytical mathe- 
matical treatments. The study of such complex flow fields requires the development 
of computational schemes which can obtain approximate numerical solutions to the 

unsteady equations governing the flow. 

Up until recently, most numerical efforts to solve these kinds of unsteady flows 
concentrated on potential methods. The early work dealt with solving the unsteady 
small disturbance potential equation as a way of obtaining the unsteady flow around 
oscillating airfoils or cascades. Later work was directed toward solving the linearized 
unsteady potential equation and the unsteady full potential equation. References 
14 through 18 represent some of the numerical work that has been done using an 

unsteady potential formulation. 

Potential methods have proven to work weU for oscillating airfoil problems, 
but unfortunately they cannot adequately account for the vortical part of the flow. 
Previous potential formulations which have included the effects of the upstream 
vorticity have invoked the linear thin airfoil approximation and assumed that the 
imposed vortical gust is convected without distortion by the nonuniform mean flow. 
This was the approach used by McCroskey and Goorjian 19 and McCroskey 20 . How- 
ever, as shown by Goldstein and Atassi 7 and Atassi 8 , the assumption that the gust 
is convected without distortion is not justified and is a poor approximation for flows 
with a spatially varying mean flow. This is especially true for turbomachinery and 
propeller flow fields where the blades are heavily loaded and there are strong mean 

flow gradients. 

In the past few years, computational efforts in unsteady aerodynamics have 
concentrated on the so-called primitive variable methods, wherein a system of non- 
linear equations such as the unsteady Euler or Navier-Stokes equations are solved 
in time along with certain specified boundary conditions. References 21 through 
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25 represent some of the recent work in this area. Unlike the potential methods, 
the primitive variable methods are equally well-suited to both oscillating airfoil 
problems and flows with convected upstream vorticity. 

The major disadvantage of the primitive variable approach is that it requires 
solving a system of unsteady, nonlinear, partial differential equations. While much 
progress has been made in this area, there are still a number of difficulties associated 
with this approach. Among the difficulties are the need for large computer mem- 
ory and lengthy solution times. Sophisticated algorithms are required to solve the 
equations, and special grid generation capabilities may also be needed. In addition, 
uncertainties about far field boundary conditions for time dependent calculations 
governed by nonlinear equations leaves some question about the accuracy of the 
solutions. More progress is needed in the development of efficient, accurate algo- 
rithms and accurate far field boundary conditions for unsteady, nonlinear problems 
to make the primitive variable approach more generally useful. 

In many unsteady flows, and in particular many turbomachinery flows, the 
upstream disturbances are small compared to the mean flow, and it is possible to 
use the linearized approach of the rapid distortion theory to study these “weakly 
rotational flows. 4 Goldstein 26 recently proposed a general theory for this approach 
which leads to a formulation which is valid for compressible, three-dimensional 
vortical flows around arbitrary bodies. Considering the problem of periodic vortical 
and entropic disturbances imposed upstream on a mean potential flow field, he 
showed that the unsteady velocity due to the imposed upstream disturbances can 
be decomposed into the sum of a known vortical component, a known entropic 
component, and an unknown irrotational (potential) component. The expressions of 
the vortical and entropic components are known functions of the upstream unsteady 
disturbances and the mean flow Lagrangian coordinates. The unsteady potential 
satisfies a nonconst ant -coefficient, inhomogeneous, convective wave equation. 

Goldstein’s decomposition of the unsteady velocity greatly simplifies the math- 
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ematical treatments of potential mean flows with imposed upstream vortical and 
entropic disturbances. Whereas a previous formulation by Hunt 27 led to a sys- 
tem of three Poisson equations for the vortical velocity and a Laplace equation for 
the potential part of the unsteady velocity, Goldstein’s formulation leads to a sin- 
gle convective wave equation for the potential part of the unsteady velocity, and 
the vortical part is obtained from its known functional expression of the upstream 
disturbances and mean flow quantities. 

For most aerodynamic flows of practical interest, there will be a stagnation 
point near the airfoil leading edge, and at this point the mean flow Lagrangian 
coordinates will become singular and will remain so along the airfoil surface and 
its wake. Since Goldstein’s vortical velocity is proportional to the gradients of 
the Lagrangian coordinates, it will also be singular along these surfaces. In order 
to ensure that the normal velocity vanishes at the airfoil surface, it will then be 
necessary for the potential part of the unsteady velocity to cancel the singular 
behavior of the vortical velocity, and therefore the potential must satisfy a singular 
airfoil boundary condition. Because of these features, it is not possible to use 
Goldstein’s formulation directly for numerical computations for aerodynamic flows 
with a stagnation point. 

Atassi and Grzedzinski 28 have shown that it is possible to modify Goldstein s 
splitting of the unsteady velocity in such a way as to remove the singular and 
indeterminate character of the resulting airfoil and wake boundary condition for the 
unsteady potential. In their decomposition of the unsteady velocity, the vortical 
component is a known function of the upstream disturbance conditions and the 
Lagrangian coordinates of the mean potential flow, and has the remarkable property 
that its normal and streamwise velocity components vanish at the surface of the 
airfoil and in its wake. The unsteady potential satisfies Goldstein’s convective wave 
equation with a modified source term. 

For incompressible flows the governing convective wave equation reduces to a 
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Poisson equation. Atassi and Grzedzinski 29 used a Green’s function formulation to 
study incompressible, vortical flows about isolated airfoils with arbitrary thickness, 
camber, and angle of attack. Their analysis assumed three-dimensional, periodic 
vortical disturbances imposed upstream of a two-dimensional airfoil of infinite span. 
Their results showed that the airfoil geometry and the upstream three-dimensional 
conditions of the vortical disturbances has a strong influence on the airfoil unsteady 
response. 

For compressible flows, the governing convective wave equation of Goldstein as 
modified by Atassi and Grzedzinski has variable coefficients which are functions of 
the mean flow quantities. In this case it is not possible to formulate the problem 
in terms of an integral equation or to derive an analytical solution. Due to the 
effects of compressiblity, therefore, it is necessary to use a full numerical approach 
to obtain approximate solutions to the resulting boundary value problem. 

The major purpose of the present work is to present a frequency domain nu- 
merical scheme which has been developed to implement this linearized approach in 
order to obtain solutions to periodic vortical flows around lifting airfoils of arbitrary 
thickness, camber, and angle of attack in compressible, subsonic flow. Previous 
papers 30 ’ 31 have presented details of the present approach for the special case of 
nonlifting airfoils. 

The linearized unsteady aerodynamic analysis and solution method that are 
presented represent an alternative to the potential and primitive variable methods 
for the solution of unsteady aerodynamic flow problems. The present method offers 
the computational efficiency of potential methods, but at the same time accounts 
for the convection and distortion of the upstream vorticity by the nonuniform mean 
flow. Our approach is therefore equally well-suited to vortical flow problems as 
well as oscillating airfoil problems. In addition, since our linearization is about the 
nonuniform mean flow, the full nonlinear effects of the mean flow are accounted for. 
Only the unsteady part of the flow is linearized. These features, coupled with the 
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inherent efficiency of the linearized approach, make the preaent approach an ideal 

solution method for unsteady aerodynamic flow fields. 

In addition to presenting the extension of our numerical scheme to lifting air- 
foils, another purpose of the present dissertation is to present numerical solutions 
for a wide range of flow conditions and airfoil geometries in order to assess the 
relative effects of Mach number, airfoil thickness, camber and angle of attack, and 
upstream disturbance conditions on the airfoil unsteady response (unsteady lift and 
moment). Our analysis will consider imposed upstream three-dimensional, vortical 
disturbances around two-dimensional airfoils of infinite span. We will assume that 
there are no incident acoustic waves or upstream entropy disturbances. 

Although we will not directly concern ourselves with the acoustic radiation 
in the present dissertation, the method that we present can also give the far field 
unsteady pressure from which the acoustic radiation can be calculated. This avoids 
the use of Lighthill’s analogy and the tedious calculation which may result from 
the refraction and scattering of the radiated acoustic waves by the spatially varying 
flow field. In a future paper, the present approach will be used to calculate the far 
field acoustic radiation for arbitrary lifting airfoils in subsonic flows. Dusey 32 , and 
Atassi and Dusey 33 , have carried out a similar calculation for flat plate airfoils. 

In Chapter II we present the general aerodynamic theory and formulate the 
boundary value problem for the linearized treatment of potential mean flows with 
imposed upstream, periodic vortical disturbances. In Chapter III we provide the 
details of the general frequency domain numerical scheme which has been developed 
for lifting airfoils in subsonic flow. Finally, we present and discuss numerical results 
for a wide variety of flow conditions in Chapter IV. 
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CHAPTER II 


AERODYNAMIC THEORY 


2.1 Linearized Euler Equations 


Consider an inviscid, compressible flow past an airfoil placed at nonzero in- 
cidence to a stream with uniform upstream velocity in the a?i direction. We 
shall assume in the present discussion that there are no shocks in the flow. Now 
if we also assume the fluid to be an ideal, non-heat conducting gas with constant 
specific heats, then the governing continuity, momentum, and entropy conservation 


equations can be written 


Dp -» 


DU - 
= ~ Vp 


Ds 

Dt 


= 0 


( 2 . 1 ) 

( 2 . 2 ) 

(2.3) 


where ^ is the material derivative, t denotes time, and p, U, p and s denote the 
fluid density, velocity, pressure, and entropy, respectively. In addition, by ideal gas, 
we mean that the pressure p, density p, and temperature T are related through a 
gas constant R by the relation p = pRT. 

In the absence of upstream flow disturbances, and assuming the airfoil to be 
rigid, there will be a steady flow U 0 (x ) about the airfoil such that U 0 (x) -► UrJ 
35 x i —* —oo, where i is a unit vector in the x\ direction. Let us suppose that far 
upstream a small amplitude, unsteady velocity disturbance, which we will denote by 
Uqo, is imposed on the flow. Then since Itfool is small compared to we assume 
that there will be small, unsteady perturbations in the physical properties of the 
otherwise steady flow. It is therefore reasonable to linearize the governing equations 
(2.1) - (2.3) about the mean flow state and to introduce perturbation quantities as 
follows: 
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€(x,i) = €o(x) + u(x,i) (2.4) 

pix,t)=po(S) + pV,t) (2-5) 

p(2,t) = Po(Z) + p'frt) ( 2 - 6 ) 

s(x,t) = so + s'(x,t) (2-7) 


where 0 subscripts denote the steady mean flow quantities, the entropy s 0 is con- 
stant, and u, p\ p\ and s' are the unsteady perturbation velocity, pressure, density 

and entropy, respectively. 

Substituting relations (2.4) - (2.7) into equations (2.1) - (2.3) and neglecting 
products of small quantities, one obtains the linearized continuity, momentum, and 
entropy conservation equations 


Dop ' + p'V • + V • ( p 0 u ) = 0 

i-/ t 

(2.8) 

w ( + 3 W„) + l>'V« ■ Vt? 0 = -Vp' 

!✓ L 

(2.9) 

y = 0 , 

(2.10) 


w here = §- t + U 0 • V is the material derivative associated with the mean flow. 

Equations (2.8) - (2.10) are the governing partial differential equations for the 
unknown perturbation quantities u,p',p', and s' for the problem of small, unsteady 
disturbances to an otherwise steady flow. The upstream disturbance u 0 0 is essen- 
tially a boundary condition which is imposed on the unsteady velocity u. That is, 
we must have u as *i ->• -oo. We shall assume in the present paper that 

{Too is the only upstream disturbance imposed on the mean flow, i.e., there are no 
imposed entropy disturbances or incident acoustic waves. 

Now from previous work concerning small amplitude velocity disturbances im- 
posed on a uniform flow, it is known that the unsteady velocity can be decomposed 
into the sum of a vortical component which has zero divergence, is purely convected, 
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and is decoupled from the pressure and any thermodynamic property, and an irrota- 
tional component which is directly related to the pressure but produces no entropy 
fluctuations. 34,26 Then far upstream the total velocity U must be of the form 

& = U 00 i + u 00 (x-iU 00 t) (2.11) 

where 

Vu oo = 0. (2.12) 

To first order, there is no pressure associated with the velocity field u^,. The 
unsteady velocity tf must satisfy 

u(x,t) ► xi qq ( x ^ooO as X\ ► oo, (2.13) 

and the pressure p' obeys 

p'(x,t) — >0 as xi — > -oo. (2.14) 

2.2 Goldstein’s Approach 

Goldstein 26 proposed a general approach for the analysis of potential mean 
flows with imposed upstream vortical and entropic disturbances which greatly sim- 
plifies the mathematical treatments of such flows. Under the conditions assumed 
on the flow in previous section, the mean flow can be expressed as the gradient of 
a potential, 

U 0 (x) = (2.15) 

In the present work we consider two-dimensional mean flow, so that U 0 (x) = 

( Iff . Iff)- 

The simplest case of potential mean flow with imposed upstream disturbances 
is the thin airfoil problem where the potential flow is simply a unif orm parallel 
flow. In this case it is advantageous to split the unsteady velocity into a vortical 
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component which is solenoidal, and an irrotational component which is expressed 
as the gradient of an unsteady potential. For the general problem of nonuniform 
flows around real airfoils, however, the splitting into solenoidal and irrotational 
components does not lead to any simplification of the mathematical formulation. 

Goldstein proposed a new approach wherein the unsteady velocity is decom- 
posed into the sum of a known vortical component and an unknown potential 

*H> 

component V<^ so that 

u(x,0 = « (/) +W. (2-16) 

The vortical component is essentially a function of the upstream disturbance 
Uoo and the mean flow Lagrangian coordinates and their spatial gradients. The 
unsteady potential satisfies a nonconstant-coefficient, inhomogeneous, convective 
wave equation 

£»( 1 M) _ . fait) = iv • (A,S <,) ). (2-17) 

Dt v co 2 Dt ’ po Po 

where is the convective derivative associated with the mean flow, and co and 
po are, respectively, the mean flow speed of sound and density. <f> is related to the 
pressure by 

p' = ( 218 > 

The vortical component u ^ is in general not solenoidal, but it does become 
solenoidal far upstream where the mean flow is uniform. For flows with no upstream 
entropy disturbances, is given by 


= u„(X - TU„t) • ( 219 ) 

The components of (X — iUoot), where X = X2, X$), are essentially La- 

grangian coordinates of the mean flow fluid particles. X2 = X2 ( ® 1 , x 2 » x 3 ) and 
X 3 = X 3 (xi,X2,x 3 ) are functionally independent integrals of the equations 


dx 1 dx 2 _ dx 3 

~U\ = Th ~ lh 


(2.20a) 
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such that 


Xi — ► X 2 and X 3 — ► x 3 as Xi — ► — 00 . (2.20b) 

For two-dimensional mean flow, we may take 

X j = (2.21) 

Poo U OO 

and 

X 3 = *3, (2-22) 

where ^0 is the stream function of the mean flow and x 3 is the spatial coordinate 
in the spanwise direction. The component Xi is defined by 


X 1 = UocA, 


(2.23) 


where A is the Lighthill “drift” function 35 , which can be expressed in terms of <£ 0 
and $0 as 

f*° 1 1 

/ (jtt - 7pr)<»», (2- 24 ) 

J- 00 Uq Uoo 


A 

A = — — n + 


U 0 2 

where the integration is carried out on 'I'o = constant. The difference in A between 
two points on a streamline is the time it takes a mean flow fluid particle to traverse 
the distance between those two points. 

To complete the formulation of the problem, it is necessary to specify boundary 
conditions for the unsteady potential <f>. At the surface of the airfoil, the normal 
velocity component must vanish, so that (j> must satisfy 


(u^ + V<£) • n = 0 


(2.25) 


or 


d<f> 

dn 


— —u 


-*{/) 


n. 


(2.26) 
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where n is the unit normal at the surface of the airfoil. In the wake, the pressure is 
continuous so that <{> must obey 

f*(A « = 0, (2.27) 

where £\<f> is the jump in <f> across the vortex sheet behind the airfoil. Far upstream, 
<f> must satisfy 

— ► 0 as —* — oo. (2.28) 

2.3 Atassi and Grzedzinski’s Decomposition 
of the Unsteady Velocity 

For most flows of practical aerodynamic interest, there will be a stagnation 
point near the leading edge of the airfoil where the mean velocity Uo vanishes. At 
this point the drift function A has a logarithmic singularity, and the right hand 
side of (2.19) then has a nonintegrable reciprocal singularity. Since A is additive, 
the right hand side of (2.19) will remain undefined along the surface of the airfoil 
and in its wake. Equation (2.26) then shows that the boundary condition for <f> is 
singular along the airfoil surface. Because of these difficulties we conclude that it 
is not possible to use Goldstein’s approach directly for numerical computations of 
unsteady vortical flows around aerodynamic bodies with a stagnation point. 

Atassi and Grzedzinski 28 have shown that it is possible to modify Goldstein’s 
splitting of the unsteady velocity field in such a way as to remove the singular 
and indeterminate character of the resulting boundary condition for the unsteady 
potential at the airfoil surface. In their decomposition of the unsteady velocity, u 
is split into the sum of a known vortical component u^ R \ which has zero normal 
and streamwise velocity components on the airfoil and in the wake, and an unkown 
potential component <f > * that satisfies equation (2.17) with a modified source term, 
so that 

u(x,t) = u^ R) + Vf* (2.29) 
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where <f>* satisfies 


£o ( ±ML 

Dt cq 2 Dt 


) - - v • po^*) = -V ■ (pofi< R) ) 

Po Po 


The vortical component tP R ' is given by 


s<") = fi<'> + 


(2.30) 


(2.31) 


where <f> is a function that satisfies 


Pq^> 

Pt 


= 0 . 


(2.32) 


As shown in Reference 28, there is no pressure associated with the velocity V</>, so 
that the vortical velocity produces no pressure fluctuations. The pressure is 
determined entirely by <j>* and is given by 


-Po(x) 


D oT 

Dt ' 


(2.33) 




In order to choose a particular function <f> that cancels the singular behavior of 
along the surface of the airfoil and in its wake, the boundary condition 


(u i,) + v})-n = 0 


(2.34) 


is imposed at the airfoil surface and in the wake. Condition (2.34) should be under- 
stood as the limit as we move close to the airfoil and the wake. Details concerning 
the construction of the function <j> can be found in Reference 28. For the important 
special case of incident harmonic velocity disturbances, specific formulas for 4> are 
presented in Reference 28, and will be discussed later in the present dissertation. 

The function vS R ^ has zero normal and streamwise velocity components at the 
airfoil surface and in the wake, so that vS R ^ satisfies 


and 


vS R) • n = 0, 

(2.35) 

f = 0, 

(2.36) 
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where n and f are the unit normal and tangent vectors. The airfoil boundary 
condition for (f)* is then 

(u< H) + V<£*)-n = 0 (2.37) 

which reduces to simply 

• n = 0. (2.38) 

In addition to satisfying the governing equation (2.30) and the airfoil boundary 
condition (2.38), the unsteady potential <f>* must also satisfy appropriate boundary 
conditions in the far field, in the wake, and at the airfoil trailing edge. 

In the far field, equations (2.29) and (2.31) together with condition (2.13) imply 

that <f>* must satisfy 

V0* — ► — V<j> as x i — ► — oo. (2.39) 

In the wake of the airfoil, <{>* is not continuous but must satisfy a jump condition 
determined by the continuity of the unsteady pressure. Applying (2.33) on each side 
of the vortex sheet behind the airfoil leads to the condition 

— (A d>*) = 0 wake (2.40) 

£>r ’ 

Finally, at the trailing edge point <j>* must be continuous in the streamwise 
direction to ensure satisfaction of the Kutta condition. 

2.4 Upstream Disturbances 

In a previous paper 30 , it has been shown that the most general upstream vorti- 
cal disturbances can be represented as the sum of three-dimensional, harmonic vor- 
ticity waves in the blade frame of reference. Because the governing equation (2.30) 
is linear, we can without loss of generality consider a single Fourier component of 
the incident disturbance, and obtain the solution to more general disturbances by 
superposition. We therefore consider incident velocity disturbances of the form 

Uoo = &*<*-*'-*> (2.41) 
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where |a| -C Uoo, k is the wave number vector which specifies the direction of 
propagation of the gust, and a and k must satisfy 


a • k — 0 


(2.42) 


in order to ensure that u*, is solenoidal (satisfies the continuity equation). 

Now condition (2.13) shows that the unsteady velocity u(x,t ) must satisfy 

u(x,<) -► as xi — ► —oo. (2.43) 

Since X 2 — - and 'J'o is the stream function of a two-dimensional mean flow, 

Poo u 00 » 7 

r 

X 2 —* x 2 + - .. ln(xi 2 + x 2 2 ) -I- constant as Xj — ► — 00 , (2.44) 

2^rt/oo 

so that we do not have X 2 —* x 2 at upstream infinity, as required by condition 
(2.20b). Equations (2.29), (2.31), and (2.39) together with (2.19) then show that 

u(x,<)^ae‘ E(xi - y “*-* J+ ^ ln( * l,+I * 2)+con,u,,t - Xa) as Xl -+ -oo, (2.45) 


so that (2.43) is not satisfied. 

However, as discussed by Atassi 8 , for a real airfoil of finite span, X 2 — ► x 2 -(- 
constant. The two-dimensional approximation is only valid then for distances that 
are small compared to the airfoil span, but large compared to the chord. Thus, 
(2.44) should be considered in this limit. 

In order to avoid difficulties with upstream conditions, the imposed upstream 
disturbances should then take the general form 




(2.46) 


where X is defined by (2.21) - (2.23). Then with this definition of u 0 c , equations 
(2.19) and (2.45) show that condition (2.13) will be satisfied. 
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2.5 The Boundary Value Problem 


For convenience we now summarize the boundary value problem that has been 
developed in the first four sections of the present chapter. From this point on we 
drop the notation which was used to distinguish between the formulation of 
Atassi and Grzedzinski in Section 2.3 and the formulation of Goldstein in Section 

2 . 2 . 

We consider a potential mean flow about a two-dimensional airfoil of infinite 
span with three-dimensional, rotational velocity disturbances of the form 

^ = ae***-*-') (2.47) 

imposed upstream. The amplitude a satisfies |a| < Uoo, and X is defined by 

equations (2.21) through (2.23). 

The unsteady potential obeys the convective wave equation 

°1<1 £^.) - — V • (p 0 v<l>) = — ' V • (po^ R) )i ( 2 - 48 ) 

Dt y c 0 2 Dt 1 po Po 

where the unsteady velocity is given by 

u(x,f) = u (/l) +W (2- 49 ) 

and the unsteady pressure is determined from 

p' = < 2 - 5 °) 

In addition, <f> must satisfy the boundary conditions 


V</> ■ n — 0 

airfoil surface 

(2.51) 


wake 

(2.52a) 

A[V<^> • n] = 0 

wake 

(2.52b) 
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V4> — + — V<j> 


as x\ 


-oo, 


(2.53) 


and be continuous in the streamwise direction at the airfoil trailing edge. Note 
that for completeness we have also included boundary condition (2.52b), which 
imposes continuity of the normal velocity across the wake. For nonlifting airfoils 
this condition is automatically satisfied, since in that case ^ is an odd function with 
respect to X 2 - For lifting airfoils, however, <f> is no longer an odd function and 
condition (2.52b) must be imposed. 

To complete the mathematical formulation of the boundary value problem, the 
explicit expression of the function 4> must be given. For a complete discussion of 
the derivation of <j> for the general problem of arbitrary upstream disturbances, the 
reader should consult Reference 28. For the problem of periodic disturbances of the 
form (2.47) imposed upstream of a single obstacle which is two-dimensional (such 
as an isolated airfoil), it is shown in Reference 28 that <f> is given by 



- Q i fc 2 1 -e «* 2X2 yk.(X-lu x t) 
1 -I- ia Q Uooki k 2 


(2.54) 


where 


a — (qj, 02 )^ 3 ) and 


,dU 0 

°o = -(- 5 —) • 

On s 


(2.55) 


Here n denotes the direction of the outward unit normal, S denotes the stagnation 
point near the airfoil leading edge, and Uo = |C/o | is the magnitude of the mean 
velocity. 

With this definition of <f>, and with the upstream velocity disturbances given 
by (2.47), the vortical velocity may then be written 


u<*> = [V(a • + V4>. 


(2.56) 


This completes the linearized mathematical formulation for the general bound- 
ary value problem of unsteady vortical flow past a lifting airfoil. 
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CHAPTER III 


NUMERICAL SCHEME 


3.1 Reformulation and Nondime nsionalization 
of the Boundary Val ue Problem 

For numerical purposes it is necessary to reformulate the boundary value prob- 
lem presented in the previous section into a form more suitable for numerical com- 
putations. Of particular concern is condition (2.53). In order to facilitate the 
implementation of the far field boundary condition, it is convenient to replace <f> by 
a function whose gradient vanishes as r-too, where r is the distance from the 

airfoil center. 

To this end, we introduce the potential functions <f>i and <j> 2 , where 

<j) = (f>\ — <f>2 (3-1) 


and <f >2 is a known function which is constructed such that 

|<^2 — <^| — ► 0 as r — > 00 . 


(3.2) 


Equation (3.1) together with conditions (2.53) and (3.2) then show that the new 
potential function <f>\ will satisfy 

v<£ J -f V(j > 2 0 as r — ► 00 . (3.3) 


The problem may then be reformulated in terms of the unknown potential <f> 1 . 
To satisfy condition (3.2), the function <j > 2 must take the form 

a 2 lfci —aifc 2 l-e"*’*’ 


h ~ k^ ai+ 1 + motfoofci h 
where the vector X e satisfies 


) e 


ik + Xe-lUoot) 


(3.4) 


\X e - X| 0 as r — ♦ 00 . 


(3.5) 
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To satisfy this condition for the general problem of vortical flows around lifting 
airfoils, we define X t as follows: 


*.,i = 


*o_ 

tfoo 


-sgn(’i'o) ^ + sgn(^ 0 )tan 


i rU 0 





(3.6) 

ii 

CN 

%> 

* 

(3.7) 

il 

(3.8) 


The expression for X e> i is obtained by making a far field expansion of Lighthill’s 
drift function A in terms of $o and 'I'o- Note that the first term in the expansion is 
just yr 11 -, and that the second term arises due to the circulation around the airfoil. 
Since the second term vanishes for airfoils with zero circulation, it is clear that the 
formulation of the source term for nonlifting airfoils is much simpler than for lifting 
airfoils. 

We also point out that the first factor in brackets on the right hand side of 
equation (3.6) is discontinuous and undefined at the points on the airfoil where 
= 0, = 0+ and $o = 0, 'I'o = 0~. The second factor in brackets is not part of 

the expansion itself, but is included to remove the discontinuity which would occur 
due to the first factor. By including the second factor in brackets and defining 
X e i = 0 at <$o = 0, 'I'o = 0, we obtain an expression for X €> \ which is everywhere 
continuous. It is important that X e< \ be continuous along the airfoil surface, for if 
it were not, the potential function (f>\ would have to satisfy a discontinuous airfoil 
boundary condition. (See condition (3.10) below.) By defining X Cj \ as in (3.6), we 
obtain an expression for X t which is everywhere continuous and also ensures that 
conditions (3.5) and (3.2) will be satisfied for both lifting and nonlifting airfoils. 

Finally, c is the airfoil chord length, and the parameter fiooi which is equal to 
y/l — Moo , where Moo is the free stream Mach number, arises due to a Gothert’s 
rule correction on the mean velocity so that the expression for X e< \ is valid for both 
compressible and incompressible flows. 
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Before presenting the reformulated boundary value problem in terms of the 
potential <f >\ , we present the nondimensionalization of the problem. We normalize 

as follows: 


Xi, X2, X3, X\ , X 2 , X^,X e> l 

by 

c 

2 

$o,r 

by 

ftfoo 

^0 

by 

2 P00U 00 

U 0 ,co 

by 

Uoo 

P 0 

by 

Poo 

P' 

by 

P00U 00 | a | 

t,A 

by 

c 

2Uoo 

LO 

by 

2V oa 

c 

k\ , ^2 ) ^.3 

by 

2 

c 


by 

m 

a 

by 

1*1 


The normalized wave number k\ = where u> and are the dimen- 

sional angular frequency and free stream velocity, respectively, is called the reduced 

frequency. 

We will assume throughout the remainder of the present chapter, unless stated 
otherwise, that all quantities are nondimensional. 

The governing equation for <j>\ is then 

_ J-v • (poWO = —v • ( P ^ R) ) 

Dt ' Co 2 Dt ’ po K Po 

+ ^(^^i)-lv.(poW 2 ) (3-9) 

Dt c 0 2 Dt po 

and the boundary conditions are 

V<i>\ • n = V </>2 • n airfoil surface (3.10) 
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- fa)) = 0 wake 

A[^(<£i — 4> 2 ) • n] = 0 wake 

V<f>i -» 0 as x\ — > — oo. 


(3.11a) 

(3.11b) 

(3.12) 


For completeness we also present the nondimensional expressions for the poten- 
tial functions <f> and <j > 2 , for the unsteady velocity and pressure, and for the upstream 
velocity disturbances. 

~ i . a 2 &i — <Ji ko 1 — e ~ tkiXt . ,r t* , 

* = + 1 < 313 > 


<f> 2 


q 2 ti-g,t 2 l-e IM| ) . 

V + l + ig<i*i *2 ' ’ l3 - 14) 


where 


Xi = A, A 2 = X e ,2 = * 0 , A 3 = X e , 3 = 


7T 


(3.15) 


= #„ - -sgn(*„)£ + «n(*„)t«i- , (- i %-)][l - e-<*S+*S>] (3,16) 

* Poo^O 


u(x,t) = u(V + V(4> 1 -<f> 2 ) 
where 


(3.17) 


it R) = [v(a • X)]^* *-*'* + V4>. 

(3.18) 

_» _ ,*\Do(4i ~ <f>i) 

p = p »(*) ■ 

(3.19) 

{Too = 

(3.20) 
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3.2 Determination of Mean Pote ntial Flow 

In order to obtain numerical solutions to equation (3.9) and its associated 
boundary conditions, one must first obtain the steady potential flow about the 
airfoil for the given flow conditions. This will in general require the use of a standard 
potential flow solver such as FL036. 36 

However, an examination of equations (3.13) through (3.18) indicate that the 
most natural choice of independent variables in which to solve equation (3.9) are 
and *o , the mean flow potential and stream functions. Since standard potential 
flow codes solve the steady problem in terms of the spatial coordinates xj and x 2 , 
there is some difficulty in obtaining the steady solution as a function of $ 0 and V 0 . 

Another difficulty arises due to the fact that the grids used by steady flow 
solvers are not suitable for the unsteady calculation. As reported in References 30 
and 31, accurate solution of equation (3.9) over a large range of flow conditions re- 
quires using grids which are determined els a function of both the reduced frequency 
Jbi and the free stream Mach number Moo • This means that in general it will be 
necessary to interpolate the solution from the steady grid onto the appropriate 
unsteady grid. 

Because of the loss of accuracy that can result from such an interpolation 
process, and also because of the need to know the mean flow as a function of $ 0 
and $ 0 , an analytical scheme that can obtain the compressible, subsonic flow about 
isolated airfoils was developed. The scheme is based on the idea that, except for a 
small inner region surrounding the airfoil, the flow gradients are not too large. Thus 
in the large outer region extending to infinity, the mean flow is essentially governed 
by a set of linear equations. As a result, one can use Gothert’s Rule, whereby the 
compressible flow about a given airfoil can be obtained from the incompressible flow 
about a similar airfoil. 

If we let a c ,#c, and 7 c denote the angle of attack, thickness ratio, and camber 
ratio of the given airfoil in a compressible flow, then the transformed airfoil for the 


25 


incompressible flow field has angle of attack, thickness ratio, and camber ratio given 
by 

<*/ = Poo «c) 

(3.21) 


<*/ = Poo <*C 
Pi = Poo Pc 
7/ = Poo 7 c 


} 


where I subscripts denote quantities from the incompressible flow field. Using 
dimensional quantities and denoting the compressible velocity by ( Uoo + «c, vc) 
at the point (x,y) and the incompressible velocity by (U a 0 + ui,vj ) at the point 
(x/, y/), the spatial coordinates and velocity in the compressible and incompressible 
planes axe related by 


and 


X = X/ 

(3.22a) 

y/ 

y = JZ 

(3.22b) 

ui 

uc = ^ 

(3.23a) 

VI 

vc = pZ- 

(3.23b) 


It is assumed here that the free stream velocity U 0 0 is aligned with the x axis, and 
that the angle of attack, thickness ratio, and camber ratio of the airfoil are such that 
the perturbation velocities axe small compared to Uoo • The potential and stream 
functions of the compressible flow field are then related to the potential and stream 
functions of the incompressible flow field by 


$0 UooX — „ 2 UooXl) 


*0 - tfooy = - tfooy/) 

Poo 


(3.24a) 

(3.24b) 


Using (3.22), and assuming that all quantities have been nondimensionalized as in 
the previous section, equations (3.24) can be rewritten 
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$/ - M^xi = 


(3.25a) 


= /3oo*0 


(3.25b) 


The problem is then, given (*o,*o), sol™ equations (3.25) for and 

then use relations (3.22) and (3.23) to determine the spatial coordinates (*,») end 
velocity components (uc, «c) of the compressible flow field. If this can be done, then 
we have the compressible flow field determined as a function of (*„ *»)■ Note that 
this assumes that we can determine (*,,»/) and («/,»/) •> functions of (*„*/)■ 
Since there is a one to one correspondence between (ir.yr) and and be 

tween („„„) and (*„*,), this is not a difficulty. At the least, it can always be 
done numerically. For the special case of Joukowski airfoils, however, it is possible 
to express the complex potential (*;,*/) in terms of the polar coordinates (r, t), 
and then obtain („.») and (.„.,) through known functional expressions of (r, »). 

For the case of incompressible flow around a Joukowski airfoil the complex 

potential is given by 


, a 2 e io “ , Tj , ,C\ , K 

il + i */ = + — + Si 1 " ( n ’ + 


where K is an arbitrary constant and 


C' = C - Co' = 


re' 9 . 


(3.26) 


(3.27) 


Here = -e + W is a complex constant, and the parameters a,e, and c‘ depend 
on the airfoil geometry. T, is the steady circulation around the airfoil and is given 

by 


Ti = — 47rasin(a/ + n) 


(3.28) 


where /z is defined by 


sin /i = 


(3.29) 


a 
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The spatial coordinates (x/, y/) are determined from r and 9 through the Joukowski 
transformation 

d 2 

xi +i yi = ( + — (3.30) 

where the parameter d satisfies 


(e + d) 2 + e 1 = a 2 . 


(3.31) 


Finally, the velocity components (u/,v/) are given by 

C 2 lC , + ae i(2ar+ ^] 


Uf — t Vf = 


- 1 . 


(3.32) 


« + i) C' 2 

Using relations (3.26) through (3.30), equations (3.25) can be expressed in 
terms of r and 9 as 


(r H ) cos(0 — a/) — [2a sin(ar/ + }i)]9 + K 


— cos 9 — e + 


cP(r cos 9 — e) 


(r cos 9 — e) 2 + (r sin 9 + e') 2 


] = 0L*o 


(3.33a) 


(r - y)sin(^-Q'/) + 2asin(a/ + //)ln(^) = /? oo ^ 0 (3.33b) 

If for given $ 0 and 'I'o we can solve equations (3.33) for r and 9, then equations 
(3.30) and (3.32) can be used to get (x/,y/) and (u/,t;/), and equations (3.22) and 
(3.23) can be used to obtain (x, y) and (u c , v c ). Once we have obtained (u c , v c ), the 
other mean flow quantities can then be obtained from Bernoulli’s law for polytropic 
gases. 

We note that, while the system of equations (3.33) is highly nonlinear in the 
unknowns r and 9 , it can be routinely solved by two-dimensional Newton iteration. 
Once a subroutine has been developed to solve equations (3.33), the compressible 
steady flow around any Joukowski airfoil can be very efficiently obtained. In ad- 
dition, the mean flow is obtained for arbitrary <k 0 and 'i'o, so that there is no 
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restriction whatsoever on the particular grid that may be used for the unsteady 
calculation. 

The only limitation in obtaining the mean potential flow by this particular 
approach is the underlying assumption that «c and vc must be small compared to 
f/oo. This means that the method will not give a good approximation in the inner 
region and particularly near the stagnation point where the perturbation velocities 
are of the same order of magnitude as However, extensive testing of this 

particular approach and comparing with the steady potential flow solver FL036 
has shown that the region of inaccuracy is very small. Figures 1 through 4 show 
Mach number comparisons between the present approximate analytical scheme and 
FL036. The comparison is made at grid points along fixed grid lines used by FL036. 
It is seen that the agreement overall is quite good, with the exception of grid points 
on the airfoil surface that are near the stagnation point. Because of this inaccuracy , 
we use FL036 to calculate the mean flow quantities along the airfoil surface itself, 
and use the approximate analytical scheme off the airfoil except in a small region 
just upstream of the stagnation point. In this region, for airfoils that have steady 
loading, the velocities are calculated using a Taylor series expansion. For airfoils 
without steady loading, the velocities are calculated from a local analytical solution 
which is patched to the outer solution. 

3.3 Frequency Domain Formulation 

An inspection of equations (3.13), (3.14), (3.18), and (3.20) indicates that the 
time dependence of the boundary value problem presented in Section 3.1 comes 
entirely through the harmonic term e - '* 1 *. It is therefore possible to make a trans- 
formation from the time domain into the frequency domain by a simple change of 
dependent variable. By transforming the problem into the frequency domain, time 
is completely eliminated from the problem and it is possible to significantly simplify 
the mathematical formulation of the boundary value problem. 
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For the case of two-dimensional mean flow, we transform into the frequency 
domain by making the following change of dependent variable: 


fa = (p e ~ ik i t + ik » x * 


(3.34) 


By including the 1 ^ 3 X 3 term in the transformation, the harmonic dependence on the 
spanwise component X 3 is also eliminated, since all of the e'*® 1 ® terms then factor 
out from each side of the equation. This is of course possible in view of (3.20) and 
(3.15). 

Before presenting the governing equation in the new dependent variable <p, we 
introduce the linear operators £ and Co to simplify the notation. 


and 


C = £»(_* :£») _ itf.( M V) 
Dt cq 2 Dt p$ 


Co = 


w2 Dl , d 2 d 2 d 2 , 

_ 2 ( 0^2 ^iTr2 ) 


~Dtl 


dil 


dx\ 


where 


(3.35) 


(3.36) 


D 0 _ d d 
Dt 0 dt 

The governing equation then takes the form 


(3.37) 


Cfa = — V -(po^ + Cfa (3.38) 

P 0 

The operator Co is essentially the operator for the thin airfoil gust response problem. 
By writing the governing equation in the following equivalent form, the left hand 
side of the equation will exhibit the character of the thin airfoil operator in the fax 
field since £ — » Co there. 


C„<h + (C - CM , = -V . fo,ti<*>) + c<h 

Po 


(3.39) 
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We then have in terms of <p 


£o<f>i + (£ ~ Co)4>i — 


+ §£ 0 +2ik ' M2 °°W 0 +(k ' M *°° ~ klM 


r 2 rl k l M2 •, TT2 9 / M \l 


(3.40) 


+ - “*•> + + ~ ? 


+ (fl 2 -^a^+n.^aV-S^S* »2.| 

+ (^ oo -^ ^ o ) a$ 2 +(1 


Equation (3.40) may be simplified further by making the following change of 
both dependent and independent variables: 


and 


V? = ipe~ iK °* 


where 


_ MC 

Plo 


(3.41a) 


(3.41b) 


$ = $ 0 (3.42a) 

$ = /?oo*o (3.42b) 


Expressing equation (3.40) in terms of ip and the new independent variables 
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$ and 'P, one gets 


£0^1 + (£ — £0)^1 — 


g-ifcit+ifcaxSg-iK'o* 



+ 


dV , ( k \ M lo 



.dip .dip d 2 tp 
+ Mip + A 2 — + A i —+M-^ 



where A^.-.A^ are functions of ($, defined by 




PL 


k 2 M 2 

us 


d ,M 2 


ik ^m 


) 


-iK 0 [M 2 ( 


dUS 


2 iki) + U$ d ^( ) 


uS?p 0, 

Po d$ 1 * 


kJK o 
PL 


P 2 u 2 


A 2 ($,*) = M 2 ( 


dU[ 




US 


Uo_dpo 
Po d<b 


+ 2iK 0 f3 2 U 2 


A 3 ($, ¥) = -/?«, 


Po d$ 


A 4 (*,*) = p 2 00 -p 2 U 2 


A 5 (<!>,*) = l-U 2 


(3.43) 


(3.44a) 


(3.44b) 


(3.44c) 


(3.44d) 


(3.44e) 
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Note that in the far field the functions A 1 ...A 5 tend to zero so that the right 
hand side of (3.43) reduces to that of the linear thin airfoil theory. 4 

To complete the frequency domain formulation, it is necessary to present the 
right hand side of the governing equation (3.39). We will refer to the right hand 
side of (3.39) as the source term and denote it by e ~ ,klt+tkaXa S. To facilitate the 
presentation of the source term, we will write the governing equation as 

Co fa +(C- C 0 ) fa = e~ ik ''+ ikaxa S (3.45a) 

or 

Co fa+(C-C 0 )fa = e~ ik ' t+ikaxa (Si + 5 2 + 5 3 - 5 4 ) (3.45b) 

where 


3X 8 5l = Yfl . {[V(S'X)}e ik *- iklt +Vt} 

po 


-,M+i*3* 8 5 2 = V • {[V(a • X)]e ,£J?-i * lt + V«fr} 


-ik^+iksxs Q _ 

b3 DV Co 2 Dt ’ 


e -'kit+ik 3 x 3 <- 4 _ J_V • (po^fa) 

P 0 


We also introduce and ^2 where 

e -ik x t+ik 9 x a <jj> = 


— ikit+tk 3 x 3 n „ 


<f2 = fa 


(3.46a) 

(3.46b) 

(3.46c) 

(3.46d) 

(3.47a) 

(3.47b) 


The expressions for S 1 ...S 4 are then 
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(3.50a) 


+ M 2 Ul 


dV 2 

a$ 2 


where 


dy>2 

d$ 0 


= ik\ 


dX e>1 

a$ 0 


¥>2 


(3.50b) 


= 1 - — 8gn(* 0 )*o«~ ( * 2+ * S) l5 + sgn(^o) tan 

OQq 7T ^ Poo * 0 


)] 


r #30 ^0 r, _ -(^o + ^o) 

nQl + PL^o 


(3.50c) 


d 2 y? 2 
^ 2 0 


= 2^1 


d 2 X e ,i 

a$ 2 



aX e>1 

d$0 



(3.50d) 


= -£ { 2sgn(«'o)(l -2^)c" ( * S+ * S) [| + sgn(* 0 ) tan 

, M~M° -(*o+’*'o) _ 2#oo^o^o r j _ e -(*2+’ I 'o)l\ 

*5 + (*o + Pl°*lT 1 (3.50e) 
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The frequency domain governing equation is then given by 


” Ax>l a $2 + 0*2 ^ 02/ 


+ + *$£+* 


dtp 


d 2 tp d 2 x p 


d$ 


d* 


d $ 2 


d * 2 


(3.52) 


= e iKo *(S 1 +S 2 + S 3 - S 4 ) 

where A^.M are defined by equations (3.44) and S 1 ...S 4 are defined by equations 
(3.48), (3.49), (3.50), and (3.51). In the far field both the coefficients Ai„.A& and 

the source term 5 X + 5 2 + 5 3 - S 4 tend to zero 80 that the e( l uation reduces to a 
Helmholtz equation. 

We conclude this section by presenting the frequency domain formulation of 
the airfoil and wake boundary conditions given by equations (3.10), (3.11a) and 
(3.11b). In terms of the new coordinates (*, * ) and the new dependent variable V>, 
equation (3.10) for the airfoil boundary condition becomes 

dip e iKo * . Qi /floor 1 - e~* a + 02 + lapoih j^x .,1 (3.53) 

0* ~~ /?oo 7T $ 1 + 

The wake boundary condition (3.11a) becomes 

(-.Jfc, + i/ 0 2 ^)[A(V>e- K ”* - ¥>2)1 = 0, (3 54a) 

where A(xpe~ ,Ko * - y> 2 ) denotes the jump in the quantity (xPe~' K °* - <^ 2 ) across 
the vortex sheet behind the airfoil. Finally, condition (3.11b) becomes 

A [V(xpe~ iKo * - v? 2 ) ■ n] = 0. (3.54b) 
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3.4 Transformation into Computational 
Coordinates and Formulation of the 
Numerical Boundary Value Problem 

Our basic numerical approach to solving equation (3.52) is to use the method 
of finite difference approximations. By discretizing the flow field and employing 
finite differences at each grid point, a large linear system of equations is obtained 
which can be solved using a matrix solver. 

Previous experience in solving equation (3.52) for the case of flat plate and 
symmetric airfoils has shown that the independent variables ($, $) axe not suitable 
computational coordinates for the gust response problem . 30 ’ 31 There axe difficulties 
in obtaining consistently accurate results over a large range of Mach numbers and 
reduced frequencies, and also problems with the implementation of fax field bound- 
ary conditions. A transformation of the independent variables is needed which not 
only provides an adequate distribution of grid points around the airfoil in the near 
field, but also provides a distribution of grid points in the fax field which is suitable 
for acoustic wave propagation and the implementation of far field, radiation type 
boundary conditions. 

In order to satisfy these requirements, we make a transformation into the elliptic 
coordinates ( 77 , £ ) with the transformation 

$ = a* cos(7D7)cosh(7rf) (3.55a) 

# = a*sin(7T77)sinh(7r£) (3.55b) 

where a* is an arbitrary constant which will be defined later. Note that in the fax 
field the elliptic coordinates reduce essentially to cylindrical coordinates, and that 
the $ — ^ plane is mapped into a semi-infinite strip in the r) — £ plane. 

With this change of variables, the governing equation takes the following form: 
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*t + Pt + J(„ 

f + 9£ 2 9So 9™ 

+ /1, j(r,,£)v> + r 1 |^ + r a ^ + T 3 |^ + r < |^+T 5 


= e'* 0 *(Si + 5 2 + S3 - 5 4 ) J(r?,0- 

Here J(r],() is the Jacobian of the transformation ( 3 . 55 ) and is given by 

J(t7,£) = 7 r 2 a* 2 [sin 2 (7Tf7) + sinh 2 ( 7 r£)] , 
and T1...T5 are functions of (77, 0 defined by 


Ti(rj,0 = a* 7T A 2 C 3 (77,0 + a*7r^oo A 3 C 4 (t 7,0 


+ a*V# 


- 7 (^ +0 %,f) 2 


72(77, £) — — a* 7 r ^2^(77, £) + a* 7 r/?oo A$C$(r), £) 


+ 


2 2/0 v4 4 C 2 (77,0 *2 2 #3 AMrhO 

a * P°° J(n, 0 r + ■%.«)* 


„ *2 2 A 4 C 3 ( 77 , 0 2 , „* 2 _ 2*2 MiMI 

%,{) = <. * 2 j(,, 0 -- +a j(„,£) 


T 4 (,,0 = a ' 2 - 2 ^ 4 ^ + a*Vft 2 


Jfo.O 




(3.56) 


(3.57) 


(3.58a) 


(3.58b) 


(3.58c) 


(3.58d) 
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T„(,,0 = -2a x 


, q *2 2 o2 AtC , h , QCi ( ti ,() 

+ 2 ° Tfef) 


and where C\...C^ are functions of (r/, £) given by 


Ci (»?,£) = j- sinh(7r0 cosh(*r 0 {a*.0 

+ 2a* 2 7r 2 [sin 2 (7rr;)cos 2 (7rT/) - sinh 2 (7r^)cos 2 (7rr?)]| 
C 2 {i),Z) = j-sin(nr})cos(wT])^J(r),() 

+ 2a* 2 7r 2 [sinh 2 (7r£) cosh 2 (7rf) - cosh 2 ( tt £ ) sin 2 (^77)] | 


C 3 ( T7,£) = sinh(7rOcos(7rr;) 


64(77, £) = cosh(7r£)sin(7r77) 


Csfa.O = cosh(ir^)sinh(7r^)|j(»7, 0 

— 2a* 2 7r 2 [sin 2 (7r77)cosh 2 (7r£) + sin 2 (7rT/) cos 2 (7tt 7)] j 


(3.58e) 


(3.59a) 


(3.59b) 


(3.59c) 


(3.59d) 


(3.59e) 
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(3.59f) 


C 6 (r},() = j-sin(irr))cos(TCT}){j(r),() 

- 2a* 2 7r 2 [sinh 2 (7r^) cosh 2 (7r£) + sinh 2 (7rf) cos 2 ^)]} 

As before, Ai...A 5 are defined by equations (3.44). 

The airfoil boundary condition (3.53) has also been transformed through the 
change of variables (3.55). Expressing condition (3.53) in the variables (??,£)> one 
gets 


dip 


= —a* 7r sin(7rr7) 


e iK °* 

0oo 




Ql^ooT 

7T 


1 — e 




$ 


+ 


Q2 + iapoik^ i ikiX e ,i 

1 + iaofci 


( 3 . 60 ) 


We now proceed to discuss the remaining boundary conditions. First, the wake 
boundary condition (3.54a) may be integrated so that it becomes 


A^e"^ 0 *) = A<^ 2 + [A(ipe~ iKo< *) t e 


&<P2t.e. 



d* 



( 3 . 61 ) 


where the subscript t.e. denotes quantities at the airfoil trailing edge. Note that, 
in general, both ip and $ are discontinuous across the wake, so that in evaluating 
A(ipe ~ tIio *), it is necessary to take into account the jump in both ip and 4>. This 
condition is imposed for grid points on the lower side of the wake. 

On the upper side of the wake, we impose condition (3.54b), which specifies 
that the normal velocity component of the unsteady velocity is continuous across 
the wake. This may be written 




v>2) = Vo ,K °*-<p 2 ) 


( 3 . 62 ) 


where “+” and ” superscripts denote above and below the wake, respectively, 
and the derivatives in (3.62) are taken to be one-sided. 
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In order to proceed further in the development of the boundary conditions 
(3.61) and (3.62) in terms of the computational coordinates (t?,£)i it is necessary 
to first discuss in more detail the transformation (3.55). First, the constant a* is 
determined from the condition that the airfoil trailing edge point on the suction side, 
where $ = a*, ^ = 0, should map into the point rj = 0, £ = 0, and the stagnation 
point, where $ = — o*, ^ = 0 should map into the point tj 1) ( 0* Then o 

must be determined from 


a* = ]r / U 0 ds (3.63a) 

2,/t.e. 

where s denotes the arc length along the airfoil surface, t.e. denotes the airfoil 
trailing edge, and s.p. denotes the stagnation point. The steady solver FL036 is 
used to locate the stagnation point, and the integration in (3.63a) is carried out 


using trapezoidal integration. 

The suction surface of the airfoil, then, corresponds to the line segment on the 
j] axis between 0 and 1, and the pressure surface corresponds to the line segment 
on the T] axis between 1 and where < 2.0 (See Figure 5.) The value of 

r]t e .~ is determined from equations (3.55a) and (3.33a), and is given by 


m.e.- 


= 2 - - cos -1 ( — [2a(cos(a + p) — sin(a + p)(2n - p)) + K 

7T U*/?So L V 

,, 2 ( , fjacosp - e) 

00 \ cos ^ c a 2 - 2 a(e' sin/i + ecos/x) + e 2 + e ,2 ^J J ’ 


(3.63b) 

where the constant K is determined from (3.33a) and (3.63a) and is given by 


K = a* Poo - 2a 
+ M, 


cos (o + p) + p sin(a + p) 


f 2 

'oo 


cP(a cos p — e) 

a cos p — e + — — - . — * 

a 2 — 2a(e' sin p + e cos /x) + e 2 + e ■* (3.63c) 


and the parameters a, d , e, e ' , and p are as defined in Section 3.2. 

The pressure side of the wake, then, corresponds to the line segments given by 
{r} t . e .- < T) < 2.0, i = 0}, and {tj = 2.0, f > 0}. The suction side of the wake 
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corresponds to the positive ( axis. The upper boundary of the r] - ( grid, which is 
given by {£ = (maxi 0 < rj < 2}, corresponds to the far field boundary in the $ - ^ 

plane. 

Because of the discontinuity in $ across the wake grid line, the grid points on 
the upper and lower sides of the wake in the physical plane (the x x - x 2 plane) 
do not coincide. This presents some difficulty in the implementation of boundary 
conditions (3.61) and (3.62), inasmuch as these conditions both specify a relation 
that must be satisfied across the wake. The difficulty can be removed, however, by 
simply using a linear averaging of t/> at the two adjacent grid points to represent 
at an arbitrary point in between. Using such a linear averaging, then, boundary 
condition (3.61), which is imposed for wake grid points on the pressure side, becomes 


1(^*0*) 


avg 




= Ay> 2 + [A(0c' i,Co *) l . e .-A V »2 t . e .]c 


tkl f*, ' u 


d* 

0 , 


(3.64) 


The discontinuity in $ across the wake also leads to a shift in the location of 
the corresponding grid points in the physical plane on opposite sides of the wake. 
Because of this shift, the last several grid points on the pressure side in the physical 
plane extend past the last suction side wake grid point, so that the above averaging 
technique cannot be employed. (See Figure 5.) However, because these last few 
points are in the far field where the mean flow is nearly uniform, the function t/> 
behaves essentially as in the case of the thin airfoil problem, and is approximately 
an odd function of By assuming t/> to be an odd function with respect to # in 
the far field, condition (3.61) becomes 


[(V>e 


-iK 0 4>i 


' = -^Av>2 - ^ 


(3.65) 
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for extra far field wake grid points on the pressure side. 

For wake grid points on the suction side, the linear averaging technique can be 
used for all points, and condition (3.62) becomes 

[U. - Vl)] ^ (3-66) 

Expressing this in rj and £, one gets 


c -iK 0 * U 0 dx/Ti + 

L-.tfo* u ° 

sinh(7r£) dr) 1 

L sin(irrj ) d£ J 


TT+ (9ft _ 9fl x 

0 [ d$ dV ’ 


(3.67) 


when the averaged values of i/> lie on the rj axis. When the averaged values of xj) lie 
on the right hand boundary of the t) — £ grid {rj = 2.0, f > 0}, (3.66) becomes 


-iKp* Uo 

sinh(7r£) dr) 


- e 


,-iKp* Uo 

sinh(7r£) dr ) - 


avg 


Vf ( 3 * 


dp? 


) 


(3.68) 


dV d<K 

In order to complete the formulation of the boundary conditions, it is only 
necessary to specify conditions at the airfoil trailing edge and in the far field. At 
the trailing edge, there are two grid points that coincide, one corresponding to the 
suction side of the airfoil and one corresponding to the pressure side. It is therefore 
necessary to impose two conditions at the trailing edge point. At the point on the 
suction side, which corresponds to (??,£) = (0,0), the Jacobian of the coordinate 
transformation (3.55) vanishes. Since the Kutta condition requires that the velocity 
at the trailing edge be finite, we are led to the requirement that 
dxp 


dr) 


= 0 at (77,0 = (0,0). 


(3.69) 


44 



At the pressure side trailing edge point, we impose the condition that the 
unsteady pressure is continuous, 


We point out that since the Jacobian vanishes at the suction side trailing edge 
point, it is not possible to directly calculate the pressure in the (»?,0 coordinates 
at that point. However, it can be calculated using a Taylor series expansion from 
neighboring points. By using this approach, condition (3.70) can be satisfied. 

In presenting the far field boundary condition, we first comment that while 
condition (3.12) expresses the mathematical requirement that V<^ - 0 at upstream 
infinity, this condition cannot be imposed throughout the far field on a boundary 
at a finite distance from the airfoil. To implement such a condition would impose a 
reflecting boundary condition which can lead to large errors in the solution. 

To correctly model the physics of the present unsteady boundary value prob- 
lem requires that the far field boundary condition be such that it allows outgoing 
acoustic waves to leave the solution domain without being reflected back into the 
computational grid. This can be accomplished, for example, by using separation of 
variables along with a series expansion for the far field solution tf> and only accepting 
those terms in the series which represent outgoing waves. The difficulty with this 
approach, however, is that it leads to a matrix which requires pivoting and therefore 
longer solution times. In addition, since i/> is not continuous across the wake, but 
the series expansion for V> is continuous everywhere, there is an incompatibility near 
the wake which can lead to a poor solution in the far field. 

An alternative to the series expansion approach is to use a Sommerfield radia- 
tion condition on the unsteady pressure. This approach avoids both the difficulties 
associated with the series expansion method and is also easier to implement. The 
Sommerfield radiation condition for the pressure is the approach used in the present 
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work, and may be written in operator notation as 



^=0 


(3.71) 


$ = R cos 0 
= iZsin0. 


Neglecting ^ terms, this reduces to 


(3.72a) 

(3.72b) 




f l/k l^OOvO / ^3 \ol y „ 

~mv~jr ) ~ ( azn * =0 




b 

Poo 


(3.73) 


This condition is applied for all grid points such that 0 < rj < 2, £ = £ n 


3.5 Numerical Method 


In the previous section we presented the transformation into computational 
coordinates and the development of the numerical boundary value problem. The 
problem to be solved numerically consists of the governing equation (3.56), and the 
boundary conditions (3.60), (3.64), (3.65), (3.67), (3.68), (3.69), (3.70) and (3.73). 
As mentioned previously, our basic numerical approach is to use the method of 
finite difference approximations, and then to solve the resulting linear system of 
equations using a matrix solver. 

The first step in obtaining numerical solutions to equation (3.56) and its as- 
sociated boundary conditions for a given flow configuration is to calculate the 
source term S and the coefficient functions A l ...A 5 . This requires the evaluation 
of Lighthill s drift function and its first and second partial derivatives with respect 


46 



to and the evaluation of the mean flow quantities and their partial derivatives 
with respect to $ and ^ at each interior grid point. The mean flow quantities 
are obtained through the scheme presented in Section 3.2, and their derivatives are 
calculated using four point differencing. 

It should be emphasized that accurate evaluation of the source term is essential 
if accurate solutions to equation (3.56) are to be obtained, and this in turn depends 
largley upon the accurate calculation of the drift function. If the drift function 
is not calculated accurately, then the source term will not tend to zero in the far 
field and the numerical scheme will become unstable. One of the major advantages 
of using the analytical scheme outlined in Section 3.2 to obtain the mean flow is 
that it can determine the mean velocities at arbitrary ($,'!'). This means that 
for fixed \&, i.e. on a given streamline, we can determine the mean velocities for 
arbitrary <£. Since evaluation of the drift function requires the integration of the 
expression ( 777 - 773 -) with respect to $ on a fixed streamline, it is very easy to 
do the numerical integrations necessary to accurately evaluate the drift function 
A. In the actual calculations, we evaluate A at a given grid point as the sum of 
an analytically determined part and a numerically determined part. The analytical 
part comes from a far field expansion for A which is given by 


A = $ 0 - ^sgn(^o) | + sgn(# 0 )tan 


0 J 


(3.74) 


This expression can be used to accurately calculate A at some point far upstream, 
and then since A is additive, the remaining portion of the integration can be done 
numerically from the upstream location to the given grid point. The numerical 
integration is done using the trapezoid rule with variable spacing in $ to ensure 
accurate resolution near the airfoil. The first and second partial derivatives of A 
are approximated using four point differencing. 

Once S and A ] ... A 5 have been calculated, they can be stored separately and 
passed to the subroutine which sets up the matrix equation to be solved. 
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We now proceed to discuss the differencing used for the governing equation and 
boundary conditions. To represent equation (3.56) we use the standard nine-point, 
central difference computational molecule which is second order accurate in r) and £. 
(Nine points are required due to the mixed derivative term.) We assume in general 
that the spacing in each direction is nonuniform. Details of the grid spacing will 
be discussed momentarily. Each of the boundary conditions (3.60), (3.69), (3.70), 
and (3.73) are implemented using four-point, one-sided differencing which is third 
order accurate for (3.60), (3.69), and (3.70), and second order accurate for (3.73). 
Boundary conditions (3.67) and (3.68) are both implemented using three-point, 
one-sided differencing which is second order accurate. 

Obtaining a numerical solution to the finite difference equations representing 
the governing equation (3.56) and its associated boundary conditions requires solv- 
ing a large matrix equation whose size is equal to twice the number of grid points. 
There are difficulties in solving this linear system of equations because the matrix 
is not block tridiagonal and does not have a regular block structure which can be 
exploited. In addition, iterative solvers have convergence problems because the di- 
agonal dominance of the matrix changes as the parameters M 0 c , k \ , and k 3 are 
varied. 

Because of these difficulties, a general purpose sparse matrix solver was devel- 
oped which stores only the nonzero entries of the matrix and can solve an arbitrary 
sparse matrix equation using Gaussian elimination. The solver basically works by 
using an ordered list to represent the nonzero entries of each row in the matrix, 
and then inserts and deletes new entries in the rows of the matrix as multiples of 
each row are added to other rows to carry out the elimination process. The only 
requirement for the solver to work is that the matrix must be arranged such that 
it remains reasonably sparse during the elimination procedure. 

The sparse solver has both a pivoting and non-pivoting feature. However, 
as pivoting during the elimination process proved to be unnecessary, the pivoting 
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feature was not used. By not pivoting during the elimination, it was possible to 
increase the storage efficiency of the solver and thereby solve larger systems of 
equations. The increased storage efficiency was gained by using a mapping function 
to map sub-blocks of the rectangular two-dimensional arrays containing the nonzero 
entries of the matrix and their associated column numbers into singly dimensioned 
arrays which contained less unused storage. By using this technique, the storage 
efficiency of the solver was increased by about 25 percent. 

The final issue to be discussed in regard to our numerical scheme is the method 
of grid determination. As reported previously 30,31 , it is not possible to use a single 
grid and obtain accurate solutions to the gust response problem for a large range of 
reduced frequencies. Rather, the unsteady grid must be determined as a function 
of both the mean flow Mach number and the reduced frequency. 

This requirement is dictated by both the accuracy of the fax field boundary 
condition (3.73), and the need to adequately model the airfoil boundary condition 
(3.63) and the wake boundary conditions (3.64) and (3.65). The accuracy of the far 
field boundary condition depends on the reduced frequency ki and free stream Mach 
number in such a way that the parameter R, where R is the distance to 

the far field boundary, should remain at least 0(1). This shows that the location of 
the outer boundary of the grid must be determined as a function of ki and M 0 0 . In 
addition, there should be enough grid points per wavelength to accurately represent 
the airfoil and wake boundary conditions. Due to the harmonic terms containing the 
parameters Kq and ki , this shows that the 77 and £ spacing have to be determined 
as a function of k\ and M 0 0 . 

The spacing in the 77 direction is constant for 0 < 77 < 1, and then changes 
slightly but constant again for 1 < 77 < ? 7 t.e.-» and finally constant again, but 
slightly different from the two previous intervals, for 77*.*.- <77 <2. The spacing on 
0 < 77 < Tjt . e .- determines the spacing on the airfoil surface. Normally the number 
of grid points in the 77 direction varies from 40 for the low frequencies up to about 
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70 for reduced frequencies of 4. 

An optimal spacing of the grid points in the wake (the f direction) turns out 

to be more difficult to achieve than the spacing on the airfoil. Numerical studies 

of the thin airfoil gust response problem showed that the optimal choice of spacing 

was 12 uniformly spaced grid points per wavelength. For the case of the thin airfoil, 

the wake boundary condition analogous to condition (3.64) is much simpler, and it 

is possible to choose the spacing of the ( grid points such that they are uniformly 

spaced along the waves. However, for the general problem of nonuniform mean 

ik x f* jX 

flows, the waves in the wake axe distorted due to the e # * « o term, and it is 
no longer possible to determine the grid point spacing such that they are uniformly 
spaced along the waves. This does not prove to be a difficulty, however, because 
in the far field the general problem of vortical flows past a lifting airfoil reduces 
essentially to a linear problem as the mean flow tends to become uniform. So we 
may determine the wake spacing as in the case of the flat plate airfoil, and in the 
far field the grid points will be nearly uniformly spaced along the waves. 

For a flat plate airfoil, the wake boundary condition which imposes the conti- 
nuity of the pressure is given, corresponding to the transformation (3.55), by 


i^-[a* coah(ir^-)-a*] 
Ywaktjj — Yt.e.C 00 


( 3 . 75 ) 


In order to have uniformly spaced grid points along the waves in the wake, the 
argument of the exponential function should vary by equal fractional increments of 
7T. To place 12 grid points per wavelength, we are then led to the requirement that 
the location of the j th £ grid point be determined from the relation 


jr(«* cosh(Tr^) - a*] = j 


( 3 . 76 ) 


Solving for £, , we get 


+ < 3 - 77 > 
This is the method for determining the £ spacing in the far field. Near the airfoil 
the above procedure leads to a spacing which is too coarse to be used. So near the 
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airfoil we use uniform £ spacing which remains constant at some value A£, until a 
point is reached such that the determined from (3.77) satisfy (j+i - (i < 
From that point on, the spacing of the f grid points is determined from (3.77). 
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CHAPTER IV 


NUMERICAL RESULTS AND DISCUSSION 


4.1 Code Validation 


Extensive efforts were taken to validate the computer codes which were de- 
veloped to implement the numerical solution procedure which was outlined in the 
previous chapter. The validation process consisted of a combination of compar- 
ing numerical results with known analytical solutions to the classical thin airfoil 
gust response problems, comparing with the second order theory of Goldstein and 
Atassi 7 and Atassi, 8 comparing with the first order numerical results of Atassi and 
Grzedzinski, 29 and calculating solutions to limiting case problems i.e., as Mach 
number, thickness, angle of attack, or camber go to zero. 

Sample computation times for the results presented varied considerably, de- 
pending on the reduced frequency, Mach number, and airfoil loading. For thin, 
unloaded airfoils, with low reduced frequency gusts, typical solution times were 
about 20 seconds per frequency on the Cray X-MP at the NASA Lewis Research 
Center. The higher frequencies for these airfoils required on the order of 60 CPU sec- 
onds per frequency, with slightly higher solution times for the higher Mach number 
cases. For symmetric, unloaded airfoils with nonzero thickness, the solution times 
ranged from about 40 seconds for the lower frequencies up to about 150 seconds for 
the highest frequencies. Finally, for loaded airfoils, the solution times ranged from 
about 250 seconds for the low frequencies up to around 1200 seconds for the highest 
reduced frequencies. No effort was made to optimize the computational efficiency 
of the scheme, as our main purpose was to validate its accuracy. 

In the results that follow, comparisons are made for one-dimensional (trans- 
verse) gusts, two-dimensional (transverse and longitudinal) gusts, and fully three- 
dimensional gusts. See Figures 6, 7, and 8. 
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The first step in the validation process was to compare numerical results with 
known solutions to the classical thin airfoil gust response problems. In Figures 9 
and 10 we present comparisons between numerical and analytical results for the 
normalized unsteady lift for vortical flows past flat plate airfoils. The normalized 
unsteady lift, commonly called the response function, is defined by 


R L (k 1 y ^ 3 ? Mqq) — 


(4.1) 


npoocUoo\a\e twt 

where V is the unsteady lift. Figure 9.a shows a comparison between numerical 
results and the Sears solution 2 for the case of a one- dimensional (transverse) gust 


in incompressible flow, and Figures 9.b and 9.c compare numerical results and an- 
alytical results obtained from a Possio solver for a one-dimensional gust at Mach 
numbers of .5 and .8. The reduced frequency values at which the comparisons are 
made range from 0 to 4.0, and are shown below the plots. The point on the real 
axis and furthest to the right corresponds to k\ — 0, and the other points along 
the curve correspond in order to the other reduced frequency values. Figures 10. a 
through 10. c compare numerical results and analytical results from a Possio solver 
for three-dimensional gusts for Mach numbers of .1, .5, and .8. The conditions on 
the gust wave number parameters are shown below the plots. As can be readily 
seen, there is excellent agreement between the numerical and analytical results. The 
only loss of accuracy occurs when both the Mach number and reduced frequency 
become large. 

In order to assess the accuracy of the present numerical scheme for vortical 
flows around thin airfoils in which the mean flow is not uniform, we compare with 
the second order theory of Goldstein and Atassi 7 and Atassi.® The results given by 
Atassi assume a zero thickness airfoil, but account for the effects of camber and angle 
of attack on the airfoil unsteady response. In Figure 11, we compare the numerically 
computed response function with the second order theory for an incompressible flow 
with a two-dimensional gust about an airfoil with an angle of attack of two degrees 
and a camber ratio of .05. The numerically computed response function is for a 6% 
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thick Joukowski airfoil, while the second order theory does not take into account the 
airfoil thickness. The reduced frequency values used for the comparison are shown 
below the plot. As can be seen, the numerical results for the 6% thick airfoil show 
a slightly larger lift at the low frequencies, but a slightly smaller lift at the higher 
frequencies for k\ up to about 3.0. As will be seen later in the discussion on the 
effects of airfoil thickness, this effect can be attributed entirely to the thickness of 
the airfoil, so that the agreement is very good for reduced frequencies ranging from 
0 up to about 3.0. For the frequencies higher than 3.0, it is not possible to make 
any firm conclusion on the accuracy of the numerical results. 

It should be mentioned at this point that the numerical results in the present 
work assume that the origin of the imposed upstream gust is at the center of the 
airfoil. Due to the representation of the gust given by (2.47), the gust will not 
appear to be centered at the airfoil center and there will be a corresponding phase 
change in the unsteady lift and moment. For a Joukowski airfoil with an imposed 
upstream gust of the form (2.47), it can be shown that the unsteady lift and moment 
should be multiplied by the phase correction factor P c given by 

p c _ C t(*l(£-W'., r -*D + M£ln2-W'„,.-y')] ( 4 . 2a ) 

in order for the gust to be referenced with respect to the airfoil center. Here 
W 0 , r + i W Q} , is the complex constant associated with the complex potential for 
the incompressible flow around the transformed Joukowski airfoil, and x' c and y' c 
are the coordinates of the airfoil center rotated through the angle of attack a. These 
parameters are given by 


W 0 , r = ecosa — e' sin a + K 


W oi = — esina — e'cosa — — — In a 

27 r 


(4.2b) 

(4.2c) 


and 


x' c = i c cosa + y c sina 


(4 2d) 
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y'c — y c cos a — x c sina 


(4.2e) 


where 

x c = 2d - 1 (4.2f) 


Vc = 


2de' 
d + 1 


(4.2*) 


and the chord length of the airfoil is 2. The parameters e,e\ a and d are as defined 
in Section 3.2 and the constant K is given by equation (3.63c). Note that for airfoils 
with zero thickness and zero mean loading, the phase correction factor P c reduces 
to just 1. 

In the results given in Figure 11, the second order theory assumed a zero 
thickness airfoil. In order to to validate the numerical scheme for airfoils with 
thickness, we compare with the first order results of Atassi and Grzedzinski. 29 In 
Figure 12, we show comparisons for one-, two-, and three-dimensional gusts for 
incompressible flows around a 12% thick Joukowski airfoil with zero degrees angle 
of attack and zero camber ratio. The reduced frequency values for the comparison 
range from .2 to 2.5 and are shown below the plots. We limit the comparison to 
this range of k \ , since this is roughly the range of validity of the Green’s function 
approach of Atassi and Grzedzinski. The agreement between the two sets of results 
is good in general. 

The final step in the validation process was to calculate the solutions to various 
limiting case problems. The limiting case of M ^ — ► 0, i.e., the incompressible case, 
was covered above where the numerical results were compared to the Sears solution. 
We now present results for the limiting cases of airfoil thickness, angle of attack, 
and camber. 

In Figures 13 through 16, we compare numerical results for the unsteady lift 
and moment about the airfoil center of a 3% thick, symmetric Joukowski airfoil with 
that of a flat plate airfoil with zero thickness. Analogous to the response function 
Rl for the unsteady lift, we define the response function Rm for the unsteady 
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moment by 


M' 

R M (k ! , *3, Moo) = ^ 

where M' is the unsteady moment about the airfoil center. Figures 13 and 14 
present results for = .1, and Figures 15 and 16 present results for = .6. 
For both Mach numbers, it is seen that the small airfoil thickness has little effect on 
the imsteady response, except for high reduced frequencies in the two-dimensional 
gust case, where the magnitude of the unsteady lift is reduced by 15 - 20 percent. 
It would appear from these results that thickness effects become more important at 
the higher frequencies for the case of the two-dimensional gust. 

Figures 17 and 18 present comparisons between results for a 12% thick sym- 
metric Joukowski airfoil at zero degrees angle of attack and one degree angle of 
attack. All plots are for a free stream Mach number of .1. As in the case of airfoil 
thickness, the strongest effect is seen in the two-dimensional gust case. However, 
here there is a significant effect both for the low and high reduced frequencies. At 
the low frequency end, the effect is primarily a reduction in the magnitude of the 
unsteady lift, while at the high end it is primarily a change in phase of the unsteady 
lift. We also point out that, in agreement with the theoretical results of Atassi, 8 
for the transverse gust case in which the gust has only an upwash component, the 
steady loading on the airfoil has virtually no effect on the unsteady lift. 

Finally, in Figures 19 through 20, we compare results for a 12% thick Joukowski 
airfoil with no camber with a 12% thick Joukowski airfoil with a camber ratio of 
.02. The free stream Mach number for all plots is .1 and the angle of attack is 
zero degrees. The effect of the small airfoil camber here is exactly analogous to the 
effect of small angle of attack above, except that here the effect is stronger due to 
the increased steady loading on the airfoil. For the airfoil with one degree angle 
of attack, the steady lift coefficient was .12, while for the 2% cambered airfoil the 
steady lift coefficient was .27. In each case, the reduction of the quasi-steady lift 
(k\ — ► 0, k 2 — > 0) for the two-dimensional gust is directly proportional to the steady 
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loading on the airfoil with a proportionality constant of 0.26. Using the theoretical 
results of Atassi reported in Reference 8, it can be shown that for zero thickness 
airfoils in a two-dimensional gust in incompressible flow, the reduction in the quasi- 
steady lift for airfoils with small camber and angle of attack is proportional to the 
steady lift coefficient with a proportionality constant of = .23. The difference 
between the numerical and theoretical values of the proportionality constant can 
be accounted for by the fact that the theoretical result does not account for the 
thickness of the airfoil. 

Before concluding the present section, the author would like to emphasize the 
significance of the method of grid determination which was outlined in Section 3.5. 
In Figure 21 we present numerical results which demonstrate the kinds of errors that 
can occur as a result of using an inappropriate grid. The results shown in these 
figures were generated without determining the grid as a function of the reduced 
frequency. For each case shown, the same grid was used for all frequencies in the 
calculation. The grid used for each Mach number was the one normally used only 
for the highest reduced frequency. By using the grid for the highest frequency, it 
was assured that there would be sufficient grid resolution to resolve the waves for 
the lower frequencies. But as can be seen, the agreement is not nearly as good as 
when the grid is determined as a function of both the Mach number and reduced 
frequency. These results show the effect of keeping the outer boundary fixed, and 
not varying it with the reduced frequency in order to ensure that the representation 
of the far field boundary condition is sufficiently accurate. 

The results in Figure 22 show the kinds of errors that can occur when the 
grid points are not suitably spaced in the far field. The grids used for the results 
in Figure 22 used an exponentially decreasing spacing which was varied to ensure 
that there were enough grid points per wavelength to adequately model the wavelike 
structure of the solution. In addition, the location of the far field boundary was also 
varied to ensure that the far field boundary condition would be sufficiently accurate. 
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However, as the results show, there are large errors in the response function curves. 
This is due to the fact that exponential spacing is not suitable for this kind of wave 
propagation problem. 

4.2 Discussion of Numerical Results 

Our major purpose in the presentation and discussion of numerical results is 
to demonstrate the effects of airfoil thickness, angle of attack, camber, and Mach 
number on the airfoil unsteady response to imposed upstream rotational velocity 
disturbances (gusts). We will examine in detail the effects of each of these pa- 
rameters on the unsteady lift and moment of airfoils subjected to one-, two-, and 
three-dimensional gusts. First we will examine the effects of airfoil thickness, then 
the effects of mean loading through both angle of attack and camber, and finally 
the effects of Mach number. 

Before presenting the unsteady results, we show plots of the mean flow Mach 
number at the airfoil surface for the various flow configurations that were considered 
in the present work. These plots are shown in Figures 23 through 43. Plots of the 
airfoils used for these calculations are shown in Figures 44 and 45. We also present 
sample computational grids used for the unsteady calculations in Figures 46 through 
61. Since for each airfoil geometry a different grid was used for every Mach number 
and reduced frequency, these represent only a very small fraction of the unsteady 
grids used to generate the results that follow. 

4.2.1 Effects of Airfoil Thickness 

on the Unsteady Response Functions 

In Figures 62 through 82, we present results that demonstrate the effects of 
airfoil thickness on the unsteady response functions. Figures 62 through 68 show 
results for the one-dimensional gust case, Figures 69 through 75 show results for the 
two-dimensional gust case, and Figures 76 through 82 show results for the three- 
dimensional gust case. For each of the three cases, we begin by comparing results 
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for unloaded flat plate airfoils versus results for unloaded, symmetric airfoils with 
a thickness ratio of .12. Comparisons are made at Mach numbers of .1, .5, .6, and 
.7. At a Mach number of .8, there is a strong shock on the 12% thick airfoil, and 
so for this case we compare the flat plate results with results for a 6% thick airfoil. 
The comparisons for the unloaded airfoils are followed by comparisons for airfoils 
with mean loading. Here we compare results between airfoils with thickness ratios 
of .06 and .12. Each airfoil has an angle of attack of 2° and camber ratio of .05. We 
compare results at Mach numbers of .1 and .5. 

An examination of the results indicates that there is a significant difference 
in the effects of airfoil thickness on the unsteady response functions depending on 
whether the gust is one-, two-, or three-dimensional, and also depending on the 
Mach number. We will discuss the one- dimensional gust case first, then the two- 
dimensional case, and finally the three-dimensional case. 

For the one- dimensional gust, it is seen that the effects of thickness at the 
lower frequencies, i.e., for k\ ranging from 0 up to about .2, is primarily to increase 
the magnitude of the unsteady response, while the phase is only slightly altered. 
The results for the unsteady moment show that for this range of frequencies the 
effects of thickness are essentially the same irrespective of the Mach number. For 
the unsteady lift, however, the effects of airfoil thickness becomes stronger with 
increasing Mach number. At a Mach number of .1, for example, there is an increase 
in the quasi-steady lift (Jki — ► 0) from the flat plate to the 12% thick airfoil of 9.1%. 
At the .5 Mach number, however, there is a 10.9% increase, and at the .6 and .7 
Mach numbers, there is a 12.5% and 15.7% increase, respectively. As the reduced 
frequency increases above about .2, the airfoil thickness has a significant effect on 
both the magnitude and phase of the unsteady response. The plots for the unsteady 
lift show that the change in phase at the highest frequencies depends strongly upon 
the Mach number. At a Mach number of .1, the phase of the unsteady lift for the 
thick airfoil lags the corresponding phase for the flat plate, but at a Mach number 
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of .5, the phases are about the same, and for higher Mach numbers, the phase of the 
lift for the thick airfoil leads the phase of the flat plate airfoil by an amount that 
increases as the Mach number becomes higher. This is in contrast to the effect of 
thickness on the phase of the high frequencies of the unsteady moment, where the 
phase of the thick airfoil always lags the phase of the flat plate airfoil, irrespective 
of Mach number. We also note that the magnitude of the unsteady moment at the 
higher frequencies is strongly affected by the airfoil thickness as the Mach number 
increases. 

The above results are all for symmetric airfoils without mean loading. In Fig- 
ures 67 and 68, we present results demonstrating the effects of airfoil thickness on 
the unsteady response for airfoils that have steady loading. As mentioned previ- 
ously, the results in these figures axe for 6% thick versus 12% thick airfoils, with 
a camber ratio of .05 and angle of attack of 2°. Results are shown for Mach num- 
bers of .1 and .5. It is seen that the unsteady lift for the loaded airfoils in the 
one-dimensional gust case is affected by airfoil thickness in an analogous way to 
the unloaded airfoils. The general trend is for the magnitude of the unsteady lift 
of the thicker airfoil to be increased at the lower frequencies, while at the higher 
frequencies the magnitude is reduced and the phase is also altered. The results 
for the unsteady moment show little change with increased airfoil thickness for the 
airfoil configurations used in these figures. 

We now turn to discuss thickness effects for the case of the two-dimensional 
gust, i.e. a transverse and longitudinal gust. We begin our discussion with the 
comparison between the flat plate results versus the thick symmetric airfoil results 
in which there is no mean loading on the airfoil (See Figures 69 - 73.) There are 
noticeable differences in the effects of airfoil thickness between the one- and two- 
dimensional gust cases. 

At a free stream Mach number of .1, it is seen that the effect of thickness is 
to increase the magnitude of the unsteady lift at the lowest frequencies by about 
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9%, while at the highest frequencies the magnitude is reduced by nearly 40% and 
the phase has been shifted by about -90°. Remarkably, there is a range of reduced 
frequencies, roughly 0.45 < Jfci < 1.6, where the thickness has virtually no effect on 
the unsteady lift and moment. 

It is interesting to note that as the Mach number increases, the effect of airfoil 
thickness at the lowest frequencies actually decreases. The reason for this unex- 
pected result is not known at the present time. Recall that for the one-dimensional 
gust case, the effects of airfoil thickness became stronger for the low reduced fre- 
quencies as the Mach number increased. 

Figures 74 and 75 show the effects of airfoil thickness for the two-dimensional 
gust case for the two loaded airfoil configurations that were described earlier. For 
the unsteady lift, the general trend at both Mach numbers is for the magnitude 
of the unsteady response function for the thicker airfoil to be increased at the low 
frequency end and decreased at the high frequency end. The unsteady moment 
plots show almost no thickness effects. These results are essentially analogous to 
the previous results for the unloaded airfoil case, except that the previous results 
showed a stronger effect due to the fact that the comparison was between a 0% 
thick and 12% thick airfoil, rather than 6% and 12% thick airfoils. 

We conclude our discussion of airfoil thickness effects by looking at the results 
for the case of three-dimensional gusts (See Figures 76 - 82). We note first that the 
effects of thickness for three-dimensional gusts are noticeably different than for one- 
and two-dimensional gusts. There are two main differences which are apparent. 

First, for the low reduced frequencies, airfoil thickness does not lead to an 
increase in the magnitude of the unsteady response functions when the gust is 
three-dimensional. For Mach numbers up to .6, there is no significant difference in 
the response function magnitude on the low frequency end. However, for the .7 and 
.8 Mach numbers, the results show that the magnitude of both the unsteady lift and 
moment are noticeably reduced for the thick airfoil for reduced frequencies ranging 
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from 0 up to about .2, This is in direct contrast to the one- and two-dimensional 
gust cases where the magnitude of the unsteady response was in general increased 
for the thicker airfoil for this range of reduced frequencies. 

The second major difference in the effects of thickness for the three-dimensional 
gust case versus the one- and two-dimensional cases is that there is a strong effect 
on both the magnitude and phase of the unsteady response functions in the mid- 
frequency range, i.e., for k\ roughly between .3 and 1.0. The results show that 
for this range of reduced frequencies, especially for the higher Mach numbers, the 
magnitude of the unsteady lift and moment is significantly increased for the thicker 
airfoil. In the one- and two-dimensional gust cases, the effects of airfoil thickness 
were strongest at the low and high reduced frequencies, rather than for the middle 
range of frequencies. 

Finally, we note that for the high reduced frequencies at all Mach numbers 
shown, the magnitude of the unsteady lift and moment is significantly reduced for 
the thicker airfoil. This is the main respect in which the three-dimensioned gust 
results for thickness effects resemble the one- and two-dimensional gust results. 

4.2.2 Effects of Mean Airfoil Loading 

on the Unsteady Response Functions 

Figures 83 through 112 present results which show the effects of mean loading 
on the unsteady response functions. All of the plots in these figures compare re- 
sponse functions for a 12% thick airfoil with no steady loading versus a 12% thick 
airfoil with steady loading due to either angle of attack or camber or both. Com- 
parisons are made at Mach numbers of .1, .5, and .6. At a Mach number of .1, the 
loaded airfoils for which the comparison is made have the following configurations: 
o = 3°, camber = 0; a = 5°, camber = 0; a = 0®, camber = .05; a — 2°, camber = 
.05; and finally a = 0°, camber = .10. At a Mach number of .5, the loaded airfoil 
configurations are: a = 5°, camber = 0; a = 0°, camber = .05; a = 2°, camber = 
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.05; and a = 0°, camber = .10. At a .6 Mach number, the only comparison made is 
for an airfoil with a = 3° and camber = 0. The choice of the airfoil configurations 
was limited by the necessity that the mean flow be free of shocks. The configura- 
tions that were chosen enable us to look at loading effects due to angle of attack 
alone, camber alone, and angle of attack and camber together. 

It is evident from the results that the effects that steady airfoil loading has on 
the unsteady response functions vanes considerably, depending upon whether the 
gust is one-, two-, or three-dimensional. It is also evident that the effect on the 
unsteady moment differs in significant ways from the effect on the unsteady lift. 

The results for the one- dimensional gust case (Figures 83 - 92) show that there 
is not much effect on the response functions due to the airfoil loading unless the 
airfoil is substantially loaded. For the results presented, only the 10% cambered 
airfoil shows a strong effect on the unsteady lift and moment. At both the .1 and 
.5 Mach numbers, the effect on the unsteady lift for the 10% cambered airfoil is to 
produce a spike in the response function curve which occurs at about fci = 3.0. The 
results for the unsteady moment also show a spike in the response function curve 
for the .5 Mach number case. The spike here occurs at k x = 2.0, and is followed by a 
portion of the curve that is concave up rather than concave down. It is interesting to 
note that these features for both the unsteady lift and moment can be seen in their 
formative stages in the plots for the .5 Mach number airfoil with 5% camber and 
2° angle of attack (See Figure 90). We conclude that for the one- dimensional gust 
case, small amounts of mean airfoil loading have a negligible effect on the unsteady 
response functions, but as the loading increases, there eventually develops a spike 
that appears in the response function curve at the higher frequencies. The spike 
point is then followed by a portion of the curve which is markedly different from 
the corresponding curve for the unloaded airfoil case. 

In contrast to the results for the one-dimensional gust case, the two-dimensional 
gust results (Figures 93 - 102) show a strong effect due to airfoil mean loading, even 
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for only lightly loaded airfoils. The effect on the unsteady lift is seen to be much 
stronger than on the unsteady moment. 

For airfoils with small loading, the effect on the unsteady lift at the low reduced 
frequencies is primarily to shift the response function in the direction of the negative 
real axis without changing the basic shape of the curve. The amount that the curve 
is shifted is given by the change in the quasi-steady lift (Jfej -4 0, k 2 -* 0) from the 
unloaded airfoil to the loaded airfoil. 

For airfoils that are more heavily loaded, the response function for the unsteady 
lift is both shifted to the left and distorted. As the steady loading on the airfoil 
increases, the distortion of the curve becomes very strong until a point is reached 
such that the response function doubles back and crosses itself and the whole char- 
acter of the curve changes. This can be seen in the results in Figure lOl.a. Note 
that the evolution of the curves leading up to this marked change in the character 
of the unsteady lift response function can be seen in Figures 99.a and lOO.a. 

In Section 4.1, it was mentioned that for airfoils with small angle of attack and 
camber, the reduction in the quasi-steady lift due to airfoil loading is proportional 
to the steady lift coefficient C L . It is interesting to note that this result also holds 
approximately for airfoils that are more substantially loaded. Denoting the change 
in the quasi-steady lift by A C'^ and the proportionality constant by — K, we have 
approximately A C' L = —K Cl, where K = .25. This compares with a value of K = 
.26 for lightly loaded airfoils. As a general rule of thumb, then, the two-dimensional 
quasi-steady lift for an airfoil with steady loading is reduced over that of an unloaded 
airfoil of the same thickness by about 25% of the value of the steady lift coefficient. 
This is true irregardless of whether the induced loading is due to airfoil angle of 
attack or camber or both. We conclude, then, that for a two-dimensional gust, 
mean loading on the airfoil leads to a significant reduction in the magnitude of the 
unsteady lift at the lowest reduced frequencies. 

Unlike the results for the unsteady lift, the results for the unsteady moment 
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show that there is a marked difference in the effects of mean airfoil loading on 
the unsteady moment response function depending on whether the induced loading 
is due to airfoil angle of attack or camber. The results indicate that when the 
steady loading on the airfoil is induced by angle of attack, there is a reduction in 
the magnitude of the unsteady moment for the reduced frequencies up to about 
.5. The reduction in the quasi-steady moment is significant, being about 19% for 
a 12% thick airfoil at a Mach number of .1 and angle of attack of 5°, and 17% 
for the same airfoil at a .5 Mach number. However, when the steady loading is 
induced by airfoil camber, there is not much effect on the unsteady moment unless 
the airfoil is substantially loaded. For example, the plots for the airfoils with 5% 
camber and no angle of attack show almost no difference between the unsteady 
moment for the loaded and unloaded airfoils (See Figures 95.b and 99.b). Only 
when the camber ratio is .10 and the Mach number is .5 do the unsteady moment 
plots show much effect due to the airfoil loading, and here it is limited to the higher 
reduced frequencies beginning around 1.3. Our conclusion is that induced loading 
due to angle of attack significantly reduces the magnitude of the unsteady moment 
for the quasi-steady case and the lower reduced frequencies up to about .5, but 
induced loading due to airfoil camber has little effect on the unsteady moment for 
the two-dimensional gust case. 

In looking at the three-dimensional gust results (Figures 103 - 112), it is evident 
that there are both similarities and differences in the effects of mean loading on the 
unsteady lift and moment of airfoils subjected to 2-D and 3-D gusts. 

First, the effects are similar in that the unsteady lift for the three-dimensional 
gust case is also sensitive to small changes in mean airfoil loading. This is partic- 
ularly true for the quasi-steady lift and the lift at the lower reduced frequencies. 
It turns out that the reduction in the quasi-steady lift from the unloaded airfoil 
to the loaded airfoil, A C' L , is again approximately proportional to the steady lift 
coefficient Cl- The proportionality constant for the 3-D results shown in Figures 
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103 - 112 is about .10, as opposed to .25 for the previously described 2-D results. 
It should be pointed out that for all of the three-dimensional results shown in these 
figures, the value of the spanwise wave number, Jfc 3 , was fixed at .4. We can therefore 
make no conclusion as to how general this value of the proportionality constant is 
for three-dimensional gust results, or as to how the proportionality constant for the 
3-D case might be effected by a change in k 3 . 

The second main similarity between the 2-D and 3-D results is that in both 
cases the unsteady moment responds differently to airfoil steady loading depending 
on whether the loading is due to angle of attack or camber. For the 2-D case, the 
unsteady moment decreased significantly for the low reduced frequencies when the 
mean loading was due to angle of attack, but showed almost no change when the 
loading was due to airfoil camber. For the 3-D case, the low frequency unsteady 
moment again decreases in response to airfoil steady loading which is induced by 
angle of attack. However, when the mean loading is induced by airfoil camber, the 
unsteady moment for the low reduced frequencies actually increases rather than 
decreasing or showing no change at all. For example, at a Mach number of .1, the 
airfoil with 5° angle of attack and no camber shows a decrease in the quasi-steady 
moment of 12% over the same airfoil with no angle of attack. On the other hand, an 
airfoil with the same thickness ratio and Mach number that has no angle of attack 
but has a camber ratio of .05 shows a 7% increase in the quasi-steady moment. (See 
Figures 104.b and 105. b.) At a Mach number of .5, the same comparison between 
the two airfoils shows am 11% decrease in the quasi-steady moment for the airfoil 
with angle of attack but no camber, but an 8% increase for the airfoil with no angle 
of attack but 5% camber. (See Figures 108.b and 109.b.) 

So the two main similarities between the 2-D and 3-D results are that in each 
case the unsteady lift is sensitive to small changes in mean loading, and in each 
case the unsteady moment responds differently to airfoil loading induced by angle 
of attack and airfoil loading induced by camber. 
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The biggest difference between the 2-D and 3-D results is that for the 3-D 
case, mean loading does not have much effect on the results for the higher reduced 
frequencies. This is in contrast to the 2-D case, where the entire range of reduced 
frequencies showed a strong sensitivity to changes in mean airfoil loading. 

4.2.3 Effects of Airfoil Mach Number 

on the Unsteady Response Functions 

We conclude our discussion of the numerical results by looking at the effects of 
Mach number on the unsteady response functions. Figures 113 through 136 show 
comparisons of the unsteady response functions for various airfoils at different Mach 
numbers. All of the results are for 12% thick airfoils. The 1-D gust results are shown 
in Figures 113 - 120, the 2-D gust results in Figures 121 - 128, and the 3-D gust 
results in Figures 129 - 136. 

An analysis of the results would seem to indicate that there are two main effects 
on the unsteady response functions due to an increase in airfoil Mach number. 

First, the magnitude of the unsteady lift and moment are significantly increased 
for the low reduced frequencies ranging from 0 up to about .2 or .3, and then signif- 
icantly decreased for reduced frequencies above about 2.0. For reduced frequencies 
ranging from about .3 up to 2.0, the effect of the higher Mach number on the un- 
steady response functions is mainly a change in phase. These trends essentially 
hold for one-, two-, and three-dimensional gusts, particularly for airfoils without 
mean loading. For airfoils that have mean loading, the trends do not hold quite 
so strongly for the three-dimensional gust case, where there is little effect on the 
magnitude of the unsteady response for the low frequency cases. The strength of 
the Mach number effects on the two-dimensional gust results is also somewhat re- 
duced for the loaded airfoil case. For example, for an unloaded 12% thick airfoil, 
the two-dimensional quasi-steady lift increases by 12.6% when the Mach number is 
increased from .1 to .5. However, for a 12% thick airfoil with 2 ° angle of attack and 
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camber ratio of .05, the two-dimensional quasi-steady lift increases by only 7.1% 
when the Mach number is increased from .1 to .5 (See Figures 121. a and 126. a). 
The one-dimensional gust results show essentially the same trends with increasing 
Mach n um ber, irrespective of whether the airfoil has mean loading or not. 

The second major effect that an increase in airfoil Mach number has on the un- 
steady response functions is that it intensifies the distortion of the response function 
curves that occurs due to mean loading. This can be seen very clearly in Figures 
125, 126, and 127, where the response functions for the higher Mach number cases 
are much more distorted than those for the lower Mach number cases. Because the 
gradients in the mean flow are stronger at the higher Mach number, the imposed 
gust undergoes a stronger distortion with increasing Mach number, and there is a 
corresponding greater distortion of the unsteady response functions. 
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CHAPTER V 


CONCLUSION 

We conclude on the basis of the numerical results presented in the previous 
chapter that the frequency domain numerical scheme developed in the present work 
can be used to solve a wide variety of unsteady vortical flow problems for a large 
range of Mach numbers and reduced frequencies. The results presented demon- 
strated the effects of airfoil thickness, camber, angle of attack, and Mach number 
on the unsteady lift and moment of Joukowski airfoils in subsonic vortical flows. It 
was seen that the effects of these parameters on the unsteady response functions 
varied considerably, depending on both the magnitude of the reduced frequency 
and the dimensionality of the imposed upstream gust. In addition, the results have 
shown that the distortion by the nonuniform mean flow of the convected upstream 
disturbances can have a very strong effect on the unsteady airfoil response. The 
results presented have thus demonstrated the importance of having a numerical 
scheme which can handle mean flow distortion effects as well as three-dimensional 
gusts for a large range of reduced frequencies. 

Another contribution of the present work is that it presents a set of benchmarK 
solutions for realistic flow configurations on which unsteady computational fluid dy- 
namics codes can be validated. In addition, it lays the foundation for the numerical 
implementation of the linearized aerodynamic approach for the analysis of unsteady 
vortical flows around airfoils used in real applications. Although we have considered 
only Joukowski airfoil geometries in the present work, we point out that since any 
airfoil geometry can be represented using complex variables and conformal mapping 
techniques, the present method can be straightforwardly extended to other airfoils. 
The current work has therefore illuminated many of the important issues to be 
resolved in the development of such a general numerical implementation. 

We conclude, then, that the present approach represents an alternative to the 
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potential and primitive variable methods for the efficient and accurate solution of 
unsteady aerodynamic flow problems. The key features of the present frequency 
domain numerical scheme are the transformation into elliptical coordinates, the 
method of grid determination, the far field radiation condition, the formulation and 
method of evaluation of the source term, and the direct sparse matrix solver. 

Finally, it is hoped that the present theoretical study of unsteady vortical flows 
around lifting airfoils will help pave the way for the development of robust, highly- 
efficient computer codes which can be used in the engineering analysis of unsteady 
aerodynamic flow fields. 
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Figure 1. Comparison of Mach number at the airfoil surfacebetween 
FL036 and the present analytical scheme. — .6, a — U 

camber = 0, thickness = .12. 


71 


ORIGINAL PAGE IS 
OF POOR QUALITY 



MACH NUMBER 


"J 


o FL036 



STREAMH1SE COORDINATE 


Figure 2. Comparison of Mach number at the airfoil surface between 

FL036 and the present analytical scheme. Moo = .6, a = 2°, 
camber = 0, thickness = .12. 
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Comparison of Mach number along a FL036 grid line which 
wraps around the airfoil. Moo = -6, a — 0°, camber — 0, 
thickness = .12. 
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Figure 4. Comparison of Mach number along a FL036 grid line which 
wraps around the airfoil. A/oo = -6, oc — 2°, camber = 0, 
thickness = .12. 
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Figure 5. Computational and physical grids. 
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Figure 6. Airfoil in a transverse gust. 



Figure 7. Airfoil in a transverse and longitudinal gust. 



Figure 8. Airfoil in a three-dimensional gust. 




REAL LIFT 


Figure 9. Comparison between the numerically computed unsteady lift and 
analytical results for a flat plate airfoil in a transverse gust at 
(a) M = 0.1, (b) M = 0.5, and (c) M = 0.8; k x = 0.0, 0.007, 
0.027, 0.062, 0.110, 0.172, 0.248, 0.338, 0.442, 0.561, 0.694, 0.842, 
1.01, 1.18, 1.38, 1.59, 1.82, 2.07, 2.33, 2.62, 2.93, 3.26, 3.62, 4.01. 
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o Present Calculation 
a Possio Solver 


REAL LIFT 


Comparison between the numerically computed unsteady lift and 
analytical results for a flat plate airfoil in a three-dimensional gust 
at (a) M = 0.1, (b) M = 0.5, and (c) M = 0.8; fc x = 0.0, 0.007 
0.027, 0.062, 0.110, 0.172, 0.248, 0.338, 0.442, 0.561, 0.694, 0.842, 
1.01, 1.18, 1.38, 1.59, 1.82, 2.07, 2.33, 2.62^2.93, 3.26, 3.62, 4.01; 
k 3 = 0.442, |a| = 1, jj* = — j, = k 2 , a - k = 0, a 2 > 0. 





Imaginary Lift 



Real Lift 


Figure 11. Comparison between the numerically computed unsteady lift 
and the second order theory for an airfoil in a transverse 
and longitudinal gust. The second order theory does not account for 
the thickness of the airfoil. M 0 0 = .1, a = 2°, camber = .05, 
thickness ratio = .06. —ai — = -7071, k i — k-i, a$ — — 0, 

ki = 0.0, 0.01, 0.03, 0.06, 0.1, 0.2, 0.3, 0.45, 0.6, 0.8, 1.0, 1.3, 

1.6, 2.0, 2.5, 3.0, 3.5, 4.0. 
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Figure 12. Comparison between the numerically computed unsteady lift and the 
first order results of Atassi and Grzedzinski for a (a) transverse gust, 
(b) transverse and longitudinal gust with —a 3 = a 2 = .7071, k\ — k 2 , 
a 3 = k 3 = 0, and (c) a three-dimensional gust with k 3 = 0.4, |a| = 1, 
fj - = — j, ki = k 2 , a - k = 0, a 2 > 0. = .1, a = 0°, camber = 0, 

thickness ratio = .12. Results shown are the complex conjugate values 
of the unsteady lift, fcj = 0.2, 0.3, 0.45, 0.6, 0.8, 1.0, 1.3, 1.6, 2.0, 2.5. 
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Figure 13. Comparison between the unsteady lift of a flat plate airfoil and a 
3 percent thick Joukowski airfoil in a (a) transverse gust, 

(b) transverse and longitudinal gust with -a! = a 2 = .7071, *i = ■ h 
a 3 = fc 3 = 0, and (c) a three-dimensional gust with k 3 = 0.4, \a\ — 1 
£2 = _I, jfcj = k 2 ,a k = 0, a 2 > 0. Moo = -1, «* = °°^ camber = 0. 
ifci = 0.0, 0.01, 0.03, 0.06, 0.1, 0.2, 0.3, 0.45, 0.6, 0.8, 1.0, 

1.3, 1.6, 2.0, 2.5, 3.0, 3.5, 4.0. 
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Figure 14. Comparison between the unsteady moment of a flat plate airfoil and 
3 percent thick Joukowski airfoil in a (a) transverse gus > 

(b) transverse and longitudinal gust with —at — a 2 — 1 * _ 2 

= 0 and (c) a three-dimensional gust with fc 3 = 0.4, |a| - i 
= _I Jfc,’ = jfc 2 , a • k = 0, a 2 > 0. Moo = .1, « = 0°, camber = 0 

k\ = 0.0, 0.01, 0.03, 0.06, 0.1, 0.2, 0.3, 0.45, 0.6, 0.8, 1.0, 

1.3, 1.6, 2.0, 2.5, 3.0, 3.5, 4.0. 
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Figure 15. Comparison between the unsteady lift of a flat plate airfoil and a 
3 percent thick Joukowski airfoil in a (a) transverse gust, 

(b) transverse and longitudinal gust with —a j = a 2 = .7071, k 3 = 
a 3 = k 3 — 0, and (c) a three-dimensional gust with k 3 = 0.4, |a| = 
^ k 3 = k 2 , a ■ k = 0, <32 >0. Afoo = .6, a = 0°, camber = 

k\ = 0.0, 0.01, 0.03, 0.06, 0.1, 0.2, 0.3, 0.45, 0.6, 0.8, 1.0, 

1.3, 1.6, 2.0, 2.5, 3.0, 3.5, 4.0. 
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Figure 16. Comparison between the unsteady moment of a flat plate airfoil and a 
3 percent thick Joukowski airfoil in a (a) transverse gust, 

(b) transverse and longitudinal gust with —a 3 = = .7071, k j = &2, 

a 3 = k 3 = 0, and (c) a three-dimensional gust with k 3 = 0.4, |a| = 1, 

= -7, it, = jb 2 , a- k = 0, a 2 > 0. = .6, a = 0°, camber = 0. 

jti = 0.0, 0.01, 0.03, 0.06, 0.1, 0.2, 0.3, 0.45, 0.6, 0.8, 1.0, 

1.3, 1.6, 2.0, 2.5, 3.0, 3.5, 4.0. 
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Figure 17. Comparison between the unsteady lift of a 12 percent thick, Joukowski 
airfoil at 0° angle of attack and 1° angle of attack for a (a) transverse 
gust, (b) transverse and longitudinal gust with — a\ = a 2 = .7071, k\ — ^ 2 , 
a 3 = Jfc 3 = 0, and (c) a three-dimensional gust with fc 3 = 0.4, |a| = 1, 
si = _I Jki = k 2 , a ■ k = 0, a 2 > 0. Moo = -1, camber = 0, thickness 
ratio = 4 .12. fcj = 0.0, 0.01, 0.03, 0.06, 0.1, 0.2, 0.3, 0.45, 0.6, 0.8, 1.0, 

1.3, 1.6, 2.0, 2.5, 3.0, 3.5, 4.0. 



86 



( a ) 


2! 



0 o 2 


(b) £ 

<£ 


CD 

<t 




k->. 


Figure 18. Comparison between the unsteady moment of a 12 percent thick, Joukowsk 
airfoil at 0° angle of attack and 1° angle of attack for a (a^ transverse^ 
gust, (b) transverse and longitudinal gust with —a\ — a 2 — -7 J ^ * 
a 3 = k 3 = 0, and (c) a three-dimensional gust with k 3 = 0.4, |a| — 1, 

= _I Jk = Jt 2 , a • fc = 0, a 2 > 0. Moo = .1, camber = 0, thickness 
ratio = 4 12. fcj = 0.0, 0.01, 0.03, 0.06, 0.1, 0.2, 0.3, 0.45, 0.6, 0.8, 1.0, 

1.3, 1.6, 2.0, 2.5, 3.0, 3.5, 4.0. 
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Figure 19. Comparison between the unsteady lift of an uncambered Joukowski 
airfoil and an airfoil with camber ratio of .02 for a (a) transverse 
gust, (b) transverse and longitudinal gust with — aj = a 2 = .7071, fcj = k 2 , 
(13 = kz = 0, and (c) a three-dimensional gust with £3 = 0.4, |a| = 1, 

= — j, fcj = k 2 , a ■ k = 0, a 2 > 0 . A/qo = . 1, a = 0°, thickness 
ratio = .12. k t = 0.0, 0.01, 0.03, 0.06, 0.1, 0.2, 0.3, 0.45, 0.6, 0.8, 1.0, 

1.3, 1.6, 2.0, 2.5, 3.0, 3.5, 4.0. 
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Figure 20. 


Comparison between the unsteady moment of an uncambered Joukov.sl 
airfoil and an airfoil with camber ratio of .02 for a (a) transveree 
gust, (b) transverse and longitudinal gust with -aj = a 2 - ' 

a 3 = Jfc 3 = 0 and (c) a three-dimensional gust with fa - 0.4, |a| - 1, 

•2 = -1 fa = Jt 2 , a- k = 0, a 2 > 0. Moo = •!,« = 0°, thickness 

ratio = 4 .12. fa = 0.0, 0.01, 0.03, 0.06, 0.1, 0.2, 0.3, 0.45, 0.6, 0.8, 1.0, 
1.3, 1.6, 2.0, 2.5, 3.0, 3.5, 4.0. 





Figure 21. Comparison between numerical results generated on a fixed grid 
and analytical results for a flat plate airfoil in a transverse gust at 
(a) M = 0.1, (b) M = 0.5, and (c) M = 0.8; ki = 0.0, 0.007, 
0.027, 0.062, 0.110, 0.172, 0.248, 0.338, 0.442, 0.561, 0.694, 0.842, 
1.01, 1.18, 1.38, 1.59, 1.82, 2.07, 2.33, 2.62, 2.93, 3.26, 3.62, 4.01. 
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Figure 22. Numerical results for the unsteady lift of a 12 percent thick, 

symmetric Joukowski airfoil in a transverse gust. Results show grid 
dependent errors due to exponentially decreasing spacing in the 
£ direction which was used in place of condition (3.70). k\ ~ 0.1, 0 
0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, 1.0, 1.1, 1.2, 1.3, 1.4, 1.5, 1.6, 1.7, 1.8 
1.9, 2.0, 2.1, 2.2, 2.3, 2.4, 2.5, 2.6, 2.7, 2.8, 2.9, 3.0. 
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Figure 23. Mean flow Mach number at the airfoil surface for a 3 percent 
thick Joukowski airfoil with M 0 0 = .1, a = 0°, and camber 
ratio = 0. 
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Figure 25. Mean flow Mach number at the airfoil surface for a 12 percent 
thick Joukowski airfoil with M , » = .l,o = 1°, and camber 


MACH NUMBER 



Figure 26. Mean flow Mach number at the airfoil surface for a 12 percent 
thick Joukowski airfoil with Moo = .1, or = 3°, and camber 
ratio = 0. 
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Figure 27. Mean flow Mach number at the airfoil surface for a 12 percent 
thick Joukowski airfoil with Mao = .1,0 = 5°, and camber 
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Figure 98 Mean flow Mach number at the airfoil surface for a 12 percent 
thick Joukowski airfoil with Moo = •!,« = 0°* 311(1 camber 
ratio = .02. 
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Figure 29. 
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Mean flow Mach number at the airfoil surface for a 12 percent 
thick Joukowski airfoil with Moo = .1,0 = 0°, and camber 
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Figure 30. Mean flow Mach number at the airfoil surface for a 12 percent 
thick Joukowski airfoil with Moo — -1, o = 2°, and camber 
ratio = .05. 
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Figure 32. Mean flow Mach number at the airfoil surface for a 6 percent 
thick Joukowski airfoil with = .1, at = 2°, and camber 
ratio = .05. 
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Figure 33. Mean flow Mach number at the airfoil surface for a 12 percent 
thick Joukowski airfoil with = .5, a = 0°, and camber 
ratio = 0. Calculation by FL036. 
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Figure 34. Mean flow Mach number at the airfoil surface for a 12 percent 
thick Joukowski airfoil with Moo = .5, « = 5°, and camber 
ratio = 0. Calculation by FL036. 
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Figure 35. Mean flow Mach number at the airfoil surface for a 12 percent 
thick Joukowski airfoil with M 0 0 = .5, a — 0°, and camber 
ratio = .05. Calculation by FL036. 
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Figure 37. Mean flow Mach number at the airfoil surface for a 6 percent 
thick Joukowski airfoil with M 0 0 — .5, a = 2°, and camber 
ratio = .05. Calculation by FL036. 
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Figure 39. Mean flow Mach number at the airfoil surface for a 3 percent 
thick Joukowski airfoil with = .6, a = 0°, and camber 
ratio = 0. Calculation by FL036. 
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Figure 40. Mean flow Mach number at the airfoil surface for a 12 percent 
thick Joukowski airfoil with Moo = .6, a — 0°, and camber 
ratio = 0. Calculation by FL036. 




Figure 41. Mean flow Mach number at the airfoil surface for a 12 percent 
thick Joukowski airfoil with = .6, a = 3°, and camber 
ratio = 0. Calculation by FL036. 
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Figure 43. Mean flow Mach number at the airfoil surface for a 6 percent 
thick Joukowski airfoil with M 0 o = .8, a = 0°, and camber 
ratio = 0. Calculation by FL036. 
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Figure 44. Joukowski airfoil geometries considered in the present work. 

(a) thickness ratio = .03, camber = 0; (b) thickness ratio = .06, 
camber = 0; (c) thickness ratio = .06, camber = .05; 

(d) thickness ratio = .12, camber = 0. 
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Figure 45. Joukowski airfoil geometries considered in the present work. 

(a) thickness ratio = .12, camber = .02; (b) thickness ratio = .12, 
camber = .05; (c) thickness ratio = .12, camber = .10. 
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Figure 47. Computational grid in the physical plane for a 12 percent 

thick Joukowski airfoil with M ^ = .1, a = 5°, camber = 0, 
and = . 1 . 
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Figure 48. Computational grid in the transformed plane for a 12 percent 
thick Joukowski airfoil with Moo = -1 , or = 5°, camber = 0, 
and ki = 1.0. 
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Figure 49. Computational grid in the physical plane for a 12 percent 

thick Joukowski airfoil with = .1, a = 5°, camber = 0, 
and k\ = 1.0. 
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Figure 50. Computational grid in the transformed plane for a 12 percent 
thick Joukowski airfoil with Moo = .1, or = 5°, camber = 0, 
and ki = 3.0. 
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Figure 51. 
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Computational grid in the physical plane for a 12 percent 
thick Joukowski airfoil with Moo = .1, a = 5°, camber = 0, 
and k] = 3.0. 
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Figure 52. Computational grid in the transformed plane for a 12 percent 
thick Joukowski airfoil with M 0 0 = .1,0 = 2°, camber = .05, 
and k\ = .1. 
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Figure 53. 



Computational grid in the transformed plane for a 12 percent 
thick Joukowski airfoil with Moo = .1, ot = 2°, camber = .05, 
and k\ = .1. 
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Figure 54. Computational grid in the transformed plane for a 12 percent 
thick Joukowski airfoil with Moo = .5, a = 2°, camber = .05, 
and ki = 1.0. 
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X2 axis 



x\ axis 


Figure 55. Computational grid in the physical plane for a 12 percent 

thick Joukowski airfoil with Moo = .5, a = 2°, camber = .05, 
and k\ = 1.0. 
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Figure 56. 


Computational grid in the transformed plane for a 12 percent 
thick Joukowski airfoil with Moo = .5, a - 2°, camber . , 

and ki = 3.0. 
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Figure 58. 


Computational grid in the transformed plane for a 12 percent 
thick Joukowski airfoil with Moo = -6, or = 3°, camber = 0, 
and k\ = 1.0. 
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Figure 59. Computational grid in the physical plane for a 12 percent 

thick Joukowski airfoil with M 0 0 = .6, a — 3°, camber = 0, 
and k\ = 1.0. 
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Figure 60. 


Computational grid in the transformed plane for a 12 percent 
thick Joukowski airfoil with Moo = .6, a = 3°, camber - 0, 
and k\ = 3.0. 
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Figure 61. Computational grid in the physical plane for a 12 percent 

thick Joukowski airfoil with Moo = .6, a = 3°, camber = 0, 
and k\ = 3.0. 
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Fieure 62. Effect of airfoil thickness on (a) the unsteady lift and (b) the 
unsteady moment of a Joukowski airfoil m a transverse gus . 
M = .1 a = 0°, camber = 0. fcj = 0.0, 0.01, 0.03, 0.06, 
O.lT 0.2, 0.3, 0.45, 0.6, 0.8, 1.0, 1.3, 1.6, 2.0, 2.5, 3.0, 3.5, 4.0. 
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Figure 63. Effect of airfoil thickness on (a) the unsteady lift and (b) the 
unsteady moment of a Joukowski airfoil in a transverse gust. 
Mqo = .5, a = 0°, camber = 0. = 0.0, 0.01, 0.03, 0.06, 

0.1, 0.2, 0.3, 0.45, 0.6, 0.8, 1.0, 1.3, 1.6, 2.0, 2.5, 3.0, 3.5, 4.0. 
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Figure 64. Effect of airfoil gust, 

unsteady moment of a Joukowsk. arrto 

^ W>1 0.46,’ oX^ W 1.3 ** 3 0. 3-5, 4... 
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Figure 65. Effect of airfoil thickness on (a) the unsteady lift and (b) the 
unsteady moment of a Joukowski airfoil in a transverse gust. 

= .7, a = 0°, camber = 0. *, = 0.0, 0.01, 0.03, 0.06, 
0.1, 0.2, 0.3, 0.45, 0.6, 0.8, 1.0, 1.3, 1.6, 2.0, 2.5, 3.0, 3.5, 4.0. 
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Figure 67. Effect of airfoil thickness on (a) the unsteady lift and (b) the 
unsteady moment of a Joukowski airfoil in a transverse gust. 
Moo — -1, a — 2°, camber = .05. k\ = 0.0, 0.01, 0.03, 0.06, 
0.1, 0.2, 0.3, 0.45, 0.6, 0.8, 1.0, 1.3, 1.6, 2.0, 2.5, 3.0, 3.5, 4.0. 
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Figure 68. Effect of airfoil thickness on (a) the unsteady lift and (b) the 
unsteady moment of a Joukowski airfoil in a transverse gust. 
Moo = .5, a = 2°, camber = .05. fci = 0.0, 0.01, 0.03, 0.06, 
0.1, 0.2, 0.3, 0.45, 0.6, 0.8, 1.0, 1.3, 1.6, 2.0, 2.5, 3.0, 3.5, 4.0. 







Figure 69. Effect of airfoil thickness on (a) the unsteady lift and (b) the 
unsteady moment of a Joukowski airfoil in a transverse and 
longitudinal gust. M 0 0 = .1, a = 0°, camber = 0. k\ = 0.0, 
0.01, 0.03, 0.06, 0.1, 0.2, 0.3, 0.45, 0.6, 0.8, 1.0, 1.3, 1.6, 2.0, 
2.5, 3.0, 3.5, 4.0. -a x = a 2 = .7071, Jbi = k 2 , a 3 = k 3 = 0. 
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Figure 70. Effect of airfoil thickness on (a) the unsteady lift and (b) the 
unsteady moment of a Joukowski airfoil in a transverse and 
longitudinal gust. M = .5, a = 0°, camber = 0. k\ = 0.0, 
0.01, 0.03, 0.06, 0.1, 0.2, 0.3, 0.45, 0.6, 0.8, 1.0, 1.3, 1.6, 2.0, 
2.5, 3.0, 3.5, 4.0. -«i = a 2 = .7071, = k 2 , a 3 = k 3 = 0. 
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Figure 71. Effect of airfoil thickness on (a) the unsteady lift and (b) the 
unsteady moment of a Joukowski airfoil in a transverse and 
longitudinal gust. = .6, a = 0°, camber = 0. k 1 = 0.0, 

0.01, 0.03, 0.06, 0.1, 0.2, 0.3, 0.45, 0.6, 0.8, 1.0, 1.3, 1.6, 2.0, 
2.5, 3.0, 3.5, 4.0. -a x = a 2 = .7071, k x = k 2 , a 3 = k 3 = 0. 
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Figure 72. Effect of airfoil thickness on (a) the unsteady lift and (b) the 
unsteady moment of a Joukowski airfoil in a transverse and 
longitudinal gust. Moo = 7, a = 0°, camber = 0. *i - 0-0, 
0.01, 0.03, 0.06, 0.1, 0.2, 0.3, 0.45, 0.6, 0.8, 1.0, 1.3, 1.6, 2.0, 
2.5, 3.0, 3.5, 4.0. -ai = a 2 = .7071, fci = *2, a 3 - «3 - 0. 
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Figure 73. Effect of airfoil thickness on (a) the unsteady lift and (b) the 
unsteady moment of a Joukowski airfoil in a transverse and 
longitudinal gust. M 0 0 = .8, a = 0°, camber = 0. = 0.0, 

0.01, 0.03, 0.06, 0.1, 0.2, 0.3, 0.45, 0.6, 0.8, 1.0, 1.3, 1.6, 2.0, 
2.5, 3.0, 3.5, 4.0. —a\ =02 = .7071, ki = & 2 , a 3 = k 3 = 0. 



IMAGINARY MOMENT IMAGINARY LIFT 


0 o 5 0 


o 


Thickness Ratio 
Thickness Ratio 




igure 74. 


Effect of airfoil thickness on (a) the unsteady lift and (b) the 
unsteady moment of a Joukowski airfoil in a transverse and 
longitudinal gust. = .1, o = V, camber = 05 - «.< 

0.01, 0.03, 0.06, 0.1, 0.2, 0.3, 0.45, 0.6, 0.8, 1.0, 1.3, 1.6, 2 0, 
2.5, 3.0, 3.5, 4.0. -a\ = a 2 = .7071, fci = a 3 = h - 
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Figure 75. Effect of airfoil thickness on (a) the unsteady lift and (b) the 
unsteady moment of a Joukowski airfoil in a transverse and 
longitudinal gust. = .5, a = 2°, camber = .05. k x = 0.0, 
0.01, 0.03, 0.06, 0.1, 0.2, 0.3, 0.45, 0.6, 0.8, 1.0, 1.3, 1.6, 2.0, 

2.5, 3.0, 3.5, 4.0. — ai = <12 — .7071, k\ = k2, = kj — 0. 






Figure 76 Effect of airfoil thickness on (a) the unsteady lift and (b) the 
unsteady moment of a Joukowski airfoil in a three-dimensional 
gust. Moo = .1, a = 0°, camber = 0. fci = 0.0, 

0.01, 0.03, 0.06, 0.1, 0.2, 0.3, 0.45, 0.6, 0.8, 1.0, 1.3, 1.6, 2.0, 
2.5, 3.0, 3.5, 4.0. k 3 = 0.4, \a\ = 1, ^ h = h, 

a ■ k = 0, a.2 > 0. 
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Figure 77. Effect of airfoil thickness on (a) the unsteady lift and (b) the 
unsteady moment of a Joukowski airfoil in a three-dimensional 
gust. Moo = -5, o = 0°, camber = 0. = 0.0, 

0.01, 0.03, 0.06, 0.1, 0.2, 0.3, 0.45, 0.6, 0.8, 1.0, 1.3, 1.6, 2.0, 
2.5, 3.0, 3.5, 4.0. h = 0.4, \a\ = 1, = -J, *i = *2, 

a ■ k = 0 , a 2 > 0 . 
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Figure 78 Effect of airfoil thickness on (a) the unsteady lift and (b) the 
unsteady moment of a Joukowski airfoil in a three-dimensional 
gust. Moo = .6, a = 0°, camber = 0. fci - 0.0, 

0.01, 0.03, 0.06, 0.1, 0.2, 0.3, 0.45, 0.6, 0.8, 1.0, 1.3, 1.6, 2.0, 
2.5, 3.0, 3.5, 4.0. h = 0.4, |a| 
a • k = 0 , a? > 0 . 


i. ? = -!> *■ = 
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Figure 79. 




Effect of airfoil thickness on (a) the unsteady lift and (b) the 
unsteady moment of a Joukowski airfoil in a three-dimensional 
gust. Moo = .7, a — 0°, camber = 0. fci = 0.0, 

0.01, 0.03, 0.06, 0.1, 0.2, 0.3, 0.45, 0.6, 0.8, 1.0, 1.3, 1.6, 2.0, 
2.5, 3.0, 3.5, 4.0. k 3 = 0.4, |a| = 1, k x = k 2 , 

a ■ k = 0, a 2 >0. 
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Figure 80. Effect of airfoil thickness on (a) the unsteady lift and (b) the 
unsteady moment of a Joukowski airfoil in a three-dimensional 
gust. Moo = .8, ot — 0°, camber = 0. fci = 0.0, 

0.01, 0.03, 0.06, 0.1, 0.2, 0.3, 0.45, 0.6, 0.8, 1.0, 1.3, 1.6, 2.0, 
2.5, 3.0, 3.5, 4.0. h = 0-4, \a\ = 1, ff = - J. k i = fc 2, 

a k = 0 , 02 > 0 . 
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Figure 81. Effect of airfoil thickness on (a) the unsteady lift and (b) the 
unsteady moment of a Joukowski airfoil in a three-dimensional 
gust. Moo = .1, a = 2°, camber = .05. Jfc, = 0.0, 0.01, 0.03, 0.06, 
0.1, 0.2, 0.3, 0.45, 0.6, 0.8, 1.0, 1.3, 1.6, 2.0, 2.5, 3.0, 3.5, 4.0. 
k 3 = 0.4, |o| = 1, g = - j, ki = k 2 , a-k = 0, a 2 > 0. 
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Figure 82. Effect of airfoil thickness on (a) the unsteady lift and (b) the 
unsteady moment of a Joukowski airfoil in a three-dimensional 
gust. Moo = .5, o = 2°, camber = .05. k 3 = 0.0, 0.01, 0.03, 0.06, 
0.1, 0.2, 0.3, 0.45, 0.6, 0.8, 1.0, 1.3, 1.6, 2.0, 2.5, 3.0, 3.5, 4.0. 
k 3 = 0.4, |a| = 1, g- = - j, ki=k 2 ,ak = 0, a 2 > 0. 
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Figure 83. Effect of mean airfoil loading on (a) the unsteady lift and (b) the 
unsteady moment of a Joukowski airfoil in a transverse gust. 

Moo = .1, thickness ratio = .12, camber = 0. ki = 0.0, 0.01, 0.03, 
0.06, 0.1, 0.2, 0.3, 0.45, 0.6, 0.8, 1.0, 1.3, 1.6, 2.0, 2.5, 3.0, 3.5, 4.0. 
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Figure 84. Effect of mean airfoil loading on (a) the unsteady lift and (b) the 
unsteady moment of a Joukowski airfoil in a transverse gust. 

Moo = .1, thickness ratio = .12, camber = 0. = 0.0, 0.01, 0.03, 

0.06, 0.1, 0.2, 0.3, 0.45, 0.6, 0.8, 1.0, 1.3, 1.6, 2.0, 2.5, 3.0, 3.5, 4.0. 
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Effect of mem airfoil loading on (a) the unsteady lift and (b) the 
unsteady moment of a Joukowski airfoil in a transverse gust. 

Moo = .1, thickness ratio = .12, a = 0°. = 0.0, 0.01, 0.03, 0.06, 

0.1, 0.2, 0.3, 0.45, 0.6, 0.8, 1.0, 1.3, 1.6, 2.0, 2.5, 3.0, 3.5, 4.0. 
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Figure 86. 




Iffect of mean airfoil loading on (a) the unsteady lift and (b) the 
nsteady moment of a Joukowski airfoil 

= .1, thickness ratio = .12. h - 0.0, 0.01, 0. » • » 

. 1 , 0 . 2 , 0 . 3 , 0 . 45 , 0 . 6 , 0 . 8 , 1 . 0 , 1 . 3 , 1 . 6 , 2 . 0 , 2 . 5 , 3 . 0 , 3 . 5 , 4 . 0 . 
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Figure 88. Effect of mean airfoil loading on (a) the unsteady lift and (b) the 
unsteady moment of a Joukowski airfoil in a transverse gust. 

Moo = .5, thickness ratio = .12, camber = 0. k i = 0.0, 0.01, 0.03, 
0.06, 0.1, 0.2, 0.3, 0.45, 0.6, 0.8, 1.0, 1.3, 1.6, 2.0, 2.5, 3.0, 3.5, 4.0. 
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Figure 89. Effect of mean airfoil loading on (a) the unsteady lift and (b) the 
unsteady moment of a Joukowski airfoil in a transverse gust. 

Moc = .5, thickness ratio = .12, a = 0°. k\ = 0.0, 0.01, 0.03, 0.06, 
0.1, 0.2, 0.3, 0.45, 0.6, 0.8, 1.0, 1.3, 1.6, 2.0, 2.5, 3.0, 3.5, 4.0. 
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90 Effect of mean airfoil loading on (a) the unsteady lift and (b) the 
unsteady moment of a Joukowski airfo.1 m a ££££*% Q 

^^“or^l.M.3, 1 5 3.0, 2.5, 3.0, 3-5, 4.0. 
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Figure 91. Effect of mean airfoil loading on (a) the unsteady lift and (b) the 
unsteady moment of a Joukowski airfoil in a transverse gust. 

Mqo = -5, thickness ratio = .12, or = 0°. ki = 0.0, 0.01, 0.03, 0.06, 
0.1, 0.2, 0.3, 0.45, 0.6, 0.8, 1.0, 1.3, 1.6, 2.0, 2.5, 3.0, 3.5, 4.0. 
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Fie,, re 92 Effect of mean airfoil loading on (a) the unsteady lift and (b) the 

unsteady moment of a Joukowski a.rfo.1 in a transverse gus^ 

Moo = .6, thickness ratio = .12, camber = 0. *r - 0.0, 0.01, 0. , 

0 06 0.1, 0.2, 0.3, 0.45, 0.6, 0.8, 1.0, 1.3, 1.6, 2.0, 2.5, 3.0, 3.5, 4.0. 
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Figure 93. Effect of mean airfoil loading on (a) the unsteady lift and (b) the 
unsteady moment of a Joukowski airfoil in a transverse and 
longitudinal gust. = .1, thickness ratio = .12, camber = 0. 

= 0.0, 0.01, 0.03, 0.06, 0.1, 0.2, 0.3, 0.45, 0.6, 0.8, 1.0, 1.3, 1.6, 
2.0, 2.5, 3.0, 3.5, 4.0. — ctj — a? = .7071, ki — & 2 , <23 = £3 = 0. 
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Figure 94. Effect of mean airfoil loading on (a) the unsteady lift and (b) the 
unsteady moment of a Joukowski airfoil in a transverse and 
longitudinal gust. Moo = »1> thickness ratio = .12, camber — 0. 

1 fcj = 0.0, 0.01, 0.03, 0.06, 0.1, 0.2, 0.3, 0.45, 0.6, 0.8, 1.0, 1.3, 1.6, 
2.0, 2.5, 3.0, 3.5, 4.0. = a 2 = .7071, fci = k 2 , a 3 = k 3 = 0. 
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Figure 95. Effect of mean airfoil loading on (a) the unsteady lift and (b) the 
unsteady moment of a Joukowski airfoil in a transverse and 
longitudinal gust. M a 0 = .1, thickness ratio = .12, a = 0°. 

Jfcj = 0.0, 0.01, 0.03, 0.06, 0.1, 0.2, 0.3, 0.45, 0.6, 0.8, 1.0, 1.3, 1.6, 
2.0, 2.5, 3.0, 3.5, 4.0. -a x = a 2 = .7071, h = fc 2 , a 3 = k 3 = 0. 
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Figure 96. Effect of mean airfoil loading on (a) the unsteady lift and (b) the 
unsteady moment of a Joukowski airfoil in a transverse and 
longitudinal gust. M 0 0 = .1, thickness ratio = .12. 

kt = 0.0, 0.01, 0.03, 0.06, 0.1, 0.2, 0.3, 0.45, 0.6, 0.8, 1.0, 1.3, 1.6, 
2.0, 2.5, 3.0, 3.5, 4.0. -a x = a 2 = .7071, fcj = k 2 , a 3 = k 3 = 0. 
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Figure 97. Effect of mean airfoil loading on (a) the unsteady lift and (b) the 
unsteady moment of a Joukowski airfoil in a transverse and 
longitudinal gust. M 0 0 = .1, thickness ratio = .12, o = 0°. 

Jtj = 0.0, 0.01, 0.03, 0.06, 0.1, 0.2, 0.3, 0.45, 0.6, 0.8, 1.0, 1.3, 1.6, 
2.0, 2.5, 3.0, 3.5, 4.0. -aj = a 2 = .7071, = fc 2 , a 3 = k 3 = 0. 
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Figure 98. Effect of mean airfoil loading on (a) the unsteady lift and (b) the 
unsteady moment of a Joukowski airfoil in a transverse and 
longitudinal gust. = .5, thickness ratio = .12, camber = 0. 
Jbi = 0.0, 0.01, 0.03, 0.06, 0.1, 0.2, 0.3, 0.45, 0.6, 0.8, 1.0, 1.3, 1.6, 
2.0, 2.5, 3.0, 3.5, 4.0. -ai = a 2 = .7071, = fa, a 3 = fa = 0. 
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Figure 99. Effect of mean airfoil loading on (a) the unsteady lift and (b) the 
unsteady moment of a Joukowski airfoil in a transverse and 
longitudinal gust. A/qo = .5, thickness ratio = .12, a = 0°. 
ki = 0.0, 0.01, 0.03, 0.06, 0.1, 0.2, 0.3, 0.45, 0.6, 0.8, 1.0, 1.3, 1.6, 
2.0, 2.5, 3.0, 3.5, 4.0. -a x = a 2 = .7071, = fc 2 , a 3 = k 3 = 0. 
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Figure 100. 


Effect of mean airfoil loading on (a) the unsteady lift and (b) the 
unsteady moment of a Joukowski airfoil in a transverse an 
longitudinal gust. Moo = -5, thickness ratio = .12 
ki = 0.0, 0.01, 0.03, 0.06, 0.1, 0.2, 0.3, 0.45, 0.6, 0.8, . , • ^ • » 
2.0, 2.5, 3.0, 3.5, 4.0. -ai = a 2 = .7071, h = k 2 , <13 - 3 - 




Figure 101. Effect of mean airfoil loading on (a) the unsteady lift and (b) the 
unsteady moment of a Joukowski airfoil in a transverse and 
longitudinal gust. M = .5, thickness ratio = .12, a = 0°. 

*i = 0.0, 0.01, 0.03, 0.06, 0.1, 0.2, 0.3, 0.45, 0.6, 0.8, 1.0, 1.3, 1.6, 
2.0, 2.5, 3.0, 3.5, 4.0. -«i = a 2 = .7071, k x = Jfc 2 , a 3 = k 3 = 0. 
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Figure 102 Effect of mean airfoil loading on (a) the unsteady lift and (b) the 
Figure momen t of a JoukoW8ki airfoil in a transverse and 

longitudinal gust. A#*, = .6, thickness ratio = A2 camber = 0 

Jbj = 0.0, 0.01, 0.03, 0.06, 0.1, 0.2, 0.3, 0.45, 0.6 0.8, 1.0, 1.3^1.6, 
2.0, 2.5, 3.0, 3.5, 4.0. -a x = a 2 = .7071, h - k 2 , a 3 3 
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Figure 103. Effect of mean airfoil loading on (a) the unsteady lift and (b) the 
unsteady moment of a Joukowski airfoil in a three-dimensional 
gust. Moo =1, thickness ratio = .12, camber = 0. Jt x = 0.0, 
0.01, 0.03, 0.06, 0.1, 0.2, 0.3, 0.45, 0.6, 0.8, 1.0, 1.3, 1.6, 2.0, 2.5, 
3.0, 3.5, 4.0. k 3 = 0.4, |<z| = 1, ^ = -j, = k 2 , a ■ Jc = 0, 

a 2 > 0. 
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Figure 104. 


Effect of mean airfoil loading on (a) the unsteady lift and (b) the 
unsteady moment of a Joukowski airfoil in a three-dimensional 
gust. Moo = .1, thickness ratio = .12, camber = 0. k\ = 0.0, 
0.01, 0.03, 0.06, 0.1, 0.2, 0.3, 0.45, 0.6, 0.8, 1.0, 1.3, F6,^2.0, 2.5, 
3.0, 3.5, 4.0. k 3 = 0.4, |a| = 1, kl = a ' k = °’ 

09 > 0 . 
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Figure 105. Effect of mean airfoil loading on (a) the unsteady lift and (b) the 
unsteady moment of a Joukowski airfoil in a three-dimensional 
gust. M 0 o = .1, thickness ratio = .12, a = 0°. ki = 0.0, 

0.01, 0.03, 0.06, 0.1, 0.2, 0.3, 0.45, 0.6, 0.8, 1.0, 1.3, 1.6, 2.0, 2.5, 
3.0, 3.5, 4.0. Jt 3 = 0.4, |a| = 1, ^ = -J, k 1 = k 2 ,a-k = 0, 
a 2 > 0. 
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Figure 106. 


Effect of mean airfoil loading on (a) the unsteady lift and (b) the 
unsteady moment of a Joukowski airfoil in a three-dimensional 
gust. Moo = -1, thickness ratio = .12. fci = 0.0, 

0.01, 0.03, 0.06, 0.1, 0.2, 0.3, 0.45, 0.6, 0.8, 1.0, 1.3, E6,_2.0, 2.5, 
3.0, 3.5, 4.0. k 3 - 0.4, |a| = 1, ff = - *i = h, a- k = 0, 


02 > 0 . 
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Figure 107. Effect of mean airfoil loading on (a) the unsteady lift and (b) the 
unsteady moment of a Joukowski airfoil in a three-dimensional 
gust. Moo — -1) thickness ratio = .12, a = 0°. k\ = 0.0, 

0.01, 0.03, 0.06, 0.1, 0.2, 0.3, 0.45, 0.6, 0.8, 1.0, 1.3, 1.6, 2.0, 2.5, 
3.0, 3.5, 4.0. k 3 = 0.4, |d| = 1, fa. _ _I, k t = k 2 , a ■ k = 0, 
a 2 > 0. 
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Figure 108. Effect of mean airfoil loading on (a) the unsteady lift and (b) the 
unsteady moment of a Joukowski airfoil in a three-dimension 
eust. Moo = -5, thickness ratio = .12, camber = 0. - U-U 

0.01, 0.03, 0.06, 0.1, 0.2, 0.3, 0.45, 0.6, 0.8, 1.0, 1.3, 1.6 , ^ 2 5 ’ 
3.0, 3.5, 4.0. k 3 = 0.4, \a\ = 1, = “l> kl ~ k<1 ' a ' k ~ ’ 
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Figure 109. Effect of mean airfoil loading on (a) the unsteady lift and (b) the 
unsteady moment of a Joukowski airfoil in a three-dimensional 
gust. Mqo = -5, thickness ratio = .12, a = 0°. k\ = 0.0, 

0.01, 0.03, 0.06, 0.1, 0.2, 0.3, 0.45, 0.6, 0.8, 1.0, 1.3, 1.6, 2.0, 2.5, 
3.0, 3.5, 4.0. k 3 = 0.4, |a| = 1, k 1= k 2 ,a k = 0, 

a 2 > 0. 
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Figure 110. 


Effect of mean airfoil loading on (a) the unsteady lift and (b) the 
unsteady moment of a Joukowski airfoil in a three-dimensional 
?ust. Moo = -5, thickness ratio = .12. ki = 0.0, 

0.01, 0.03, 0.06, 0.1, 0.2, 0.3, 0.45, 0.6, 0.8, 1.0, 1.3, L.6, jl.0, 2.5, 
3.0, 3.5, 4.0. k 3 = 0.4, |a| = 1, ff = - *i =k 2 , a-k = 0, 

0.2 > 0 . 
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Figure 111. Effect of mean airfoil loading on (a) the unsteady lift and (b) the 
unsteady moment of ft Joukowski airfoil in a three-dimensional 
gust. Moo = -5) thickness ratio = .12, a = 0°. fcj = 0.0, 

0.01, 0.03, 0.06, 0.1, 0.2, 0.3, 0.45, 0.6, 0.8, 1.0, 1.3, 1.6, 2.0, 2.5, 
3.0, 3.5, 4.0. k 3 = 0.4, \a\ = 1, ff = k 1 = k 2 ,ak = 0, 
a 2 > 0. 
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Figure 112. Effect of mean airfoil loading on (a) the unsteady lift and (b) the 
unsteady moment of a Joukowski airfoil in a three-dimensional 
gust. Moo = -6, thickness ratio = .12, camber = 0. ki = 0.0, 
0.01, 0.03, 0.06, 0.1, 0.2, 0.3, 0.45, 0.6, 0.8, 1.0, 1.3, 1-6, _2.0, 2.5, 
3.0, 3.5, 4.0. k 3 = 0.4, \a\ = 1, = -j, = k 2 , a-k = 0, 

02 > 0 . 
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Figure 113. Effect of Mach number on (a) the unsteady lift and (b) the 
unsteady moment of a Joukowski airfoil in a transverse gust 
a = 0°, thickness ratio = .12, camber = 0. ki = 0.0, 0.01, 
0.03, 0.06, 0.1, 0.2, 0.3, 0.45, 0.6, 0.8, 1.0, 1.3, 1.6, 2.0, 2.5, 
3.0, 3.5, 4.0. 
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Figure 114. Effect of Mach number on (a) the unsteady lift and (b) the 
unsteady moment of a Joukowski airfoil in a transverse gust 
a = 0°, thickness ratio = .12, camber = 0. fci = 0.0, 0.01, 
0.03, 0.06, 0.1, 0.2, 0.3, 0.45, 0.6, 0.8, 1.0, 1.3, 1.6, 2.0, 2.5, 
3.0, 3.5, 4.0. 





Figure 115. Effect of Mach number on (a) the unsteady lift and (b) the 
unsteady moment of a Joukowski airfoil in a transverse gust 
o = 0°, thickness ratio = .12, camber = 0. fci = 0.0, 0.01, 
0.03, 0.06, 0.1, 0.2, 0.3, 0.45, 0.6, 0.8, 1.0, 1.3, 1.6, 2.0, 2.5, 
3.0, 3.5, 4.0. 
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Figure 116. Effect of Mach number on (a) the unsteady lift and (b) the 
unsteady moment of a Joukowski airfoil in a transverse gust 
o = 5°, thickness ratio = .12, camber = 0. fci = 0.0, 0.01, 
0.03, 0.06, 0.1, 0.2, 0.3, 0.45, 0.6, 0.8, 1.0, 1.3, 1.6, 2.0, 2.5, 
3.0, 3.5, 4.0. 






Figure 117. Effect of Mach number on (a) the unsteady lift and (b) the 
unsteady moment of a Joukowski airfoil in a transverse gust. 
a = 0°, thickness ratio = .12, camber = .05. k\ = 0.0, 0.01, 
0.03, 0.06, 0.1, 0.2, 0.3, 0.45, 0.6, 0.8, 1.0, 1.3, 1.6, 2.0, 2.5, 
3.0, 3.5, 4.0. 
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Figure 118. 


Effect of Mach number on (a) the unsteady lift and (b) the 
unsteady moment of a Joukowski airfoil in a transverse gust, 
a = 2°, thickness ratio = .12, camber = .05. h = 0.0, 0.01, 
0.03, 0.06, 0.1, 0.2, 0.3, 0.45, 0.6, 0.8, 1.0, 1.3, 1.6, 2.0, 2.5, 
3.0, 3.5, 4.0. 
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Figure 119. Effect of Mach number on (a) the unsteady lift and (b) the 
unsteady moment of a Joukowski airfoil in a transverse gust 
a — 0°, thickness ratio = .12, camber = .10. = 0.0, 0.0! 

0.03, 0.06, 0.1, 0.2, 0.3, 0.45, 0.6, 0.8, 1.0, 1.3, 1.6, 2.0, 2.5, 
.3.0, 3.5, 4.0. 
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Fieure 120. Effect of Mach number on (a) the unsteady lift and (b) the 
unsteady moment of a Joukowski airfoil in a transverse gust 
a = 3°, thickness ratio = .12, camber = 0. h = 0.0, 0.01, 
0.03, 0.06, 0.1, 0.2, 0.3, 0.45, 0.6, 0.8, 1.0, 1.3, 1.6, 2.0, 2.5, 
3.0, 3.5, 4.0. 
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Effect of Mach number on (a) the unsteady lift and (b) the 
unsteady moment of a Joukowski airfoil in a transverse and 
longitudinal gust, thickness ratio = .12, a = 0°, camber = 
ki = 0.0, 0.01, 0.03, 0.06, 0.1, 0.2, 0.3, 0.45, 0.6, 0.8, 1.0, 
1.3, 1.6, 2.0, 2.5, 3.0, 3.5, 4.0. -a, = a 2 = .7071, = fc 2 , 
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Figure 122 Effect of Mach number on (a) the unsteady lift and (b) the 
unsteady moment of a Joukowski airfoil in a transverse and 
longitudinal gust, thickness ratio = .12, a = 0°, camber - 0 

jfc, = 0.0, 0.01, 0.03, 0.06, 0.1, 0.2, 0.3, 0.45, 0.6, 0.8, 1.0, 

1.3 1.6, 2.0, 2.5, 3.0, 3.5, 4.0. -ai = a 2 = .7071, fci = k 2 , 





Figure 123. Effect of Mach number on (a) the unsteady lift and (b) the 
unsteady moment of a Joukowski airfoil in a transverse and 
longitudinal gust, thickness ratio = .12, a = 0°, camber = 0. 
Arj = 0.0, 0.01, 0.03, 0.06, 0.1, 0.2, 0.3, 0.45, 0.6, 0.8, 1.0, 

1.3, 1.6, 2.0, 2.5, 3.0, 3.5, 4.0. -oi = a 2 = .7071, k x = k 2 , 

«3 = = 0. 
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Figure 124. Effect of Mach number on (a) the unsteady lift and (b) the 
unsteady moment of a Joukowski airfoil in a transverse and 
longitudinal gust, thickness ratio = .12, a = 5°, camber = 0. 
jbj = 0.0, 0.01, 0.03, 0.06, 0.1, 0.2, 0.3, 0.45, 0.6, 0.8, 1.0, 

1.3, 1.6, 2.0, 2.5, 3.0, 3.5, 4.0. -m = a 2 = -7071, k i = k 2 , 
a$ = k: 





Figure 125. Effect of Mach number on (a) the unsteady lift and (b) the 
unsteady moment of a Joukowski airfoil in a transverse and 
longitudinal gust, thickness ratio = .12, a = 0°, camber = .05. 
it, = 0.0, 0.01, 0.03, 0.06, 0.1, 0.2, 0.3, 0.45, 0.6, 0.8, 1.0, 

1.3, 1.6, 2.0, 2.5, 3.0, 3.5, 4.0. -a t = a 2 = .7071, k t = k 2 , 
a 3 = k 3 = 0. 
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Figure 12G. Effect of Mach number on (a) the unsteady lift and (b) the 
unsteady moment of a Joukowski airfoil in a transverse and 
longitudinal gust, thickness ratio = .12, a = 2°, camber = .05. 
Jt, = 0.0, 0.01, 0.03, 0.06, 0.1, 0.2, 0.3, 0.45, 0.6, 0.8, 1.0, 

1.3, 1.6, 2.0, 2.5, 3.0, 3.5, 4.0. -a x = a 2 = .7071, k 1 = k 2 , 

a 3 = k 3 = 0 . 
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Effect of Mach number on (a) the unsteady lift and (b) the 
unsteady moment of a Joukowski airfoil in a transverse and 
longitudinal gust, thickness ratio = .12, a = 0°, camber = 
ki = 0.0, 0.01, 0.03, 0.06, 0.1, 0.2, 0.3, 0.45, 0.6, 0.8, 1.0, 
1.3, 1.6, 2.0, 2.5, 3.0, 3.5, 4.0. -a, = a 2 = .7071, i t, = k 2y 
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Figure 128. Effect of Mach number on (a) the unsteady lift and (b) the 
unsteady moment of a Joukowski airfoil in a transverse and 
longitudinal gust, thickness ratio = .12, a = 3°, camber = 0. 
k x = 0.0, 0.01, 0.03, 0.06, 0.1, 0.2, 0.3, 0.45, 0.6, 0.8, 1.0, 

1.3, 1.6, 2.0, 2.5, 3.0, 3.5, 4.0. -ai = a 2 = .7071, k x = k 2 , 

03 = ^3 = 0 . 
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Effect of Mach number on (a) the unsteady lift and (b) the 
unsteady moment of a Joukowski airfoil in a three-dimensional 
gust, a = 0°, thickness ratio = .12, camber = 0. k\ = 0.0, 
0.01, 0.03, 0.06, 0.1, 0.2, 0.3, 0.45, 0.6, 0.8, 1.0, 1.3, 1.6, 2.0, 2.5, 
3.0, 3.5, 4.0. k 3 = 0.4, \a\ = 1, ff = - J, k x =k 2 ,a-k = 0, 

09 > 0 . 
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Effect of Mach number on (a) the unsteady lift and (b) the 
unsteady moment of a Joukowski airfoil in a three-dimensional 
gust, a = 0°, thickness ratio = .12, camber = 0. &i = 0.0, 
0.01, 0.03, 0.06, 0.1, 0.2, 0.3, 0.45, 0.6, 0.8, 1.0, 1.3, 1.6, ^2.0, 2.5, 
3.0, 3.5, 4.0. k 3 — 0.4, |a| = 1, ^ = - j, k x = k 2 , a ■ k = 0, 
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Effect of Mach number on (a) the unsteady lift and (b) the 
unsteady moment of a Joukowski airfoil in a three-dimensional 
gust, a = 0°, thickness ratio = .12, camber = 0. k\ = 0.0, 
0.01, 0.03, 0.06, 0.1, 0.2, 0.3, 0.45, 0.6, 0.8, 1.0, 1.3, 1.6, 2.0, 2.5 
3.0, 3.5, 4.0. k 3 = 0.4, \S\ = 1, h =k 2 ,a-k = 0, 

a 2 ^ 





Figure 13 9 Effect of Mach number on (a) the unsteady lift and (b) the 

unsteady moment of a Joukowski airfoil in a three-dimensional 
gust, o = 5°, thickness ratio = .12, camber = 0. k x - 0.0, 
0.01, 0.03, 0.06, 0.1, 0.2, 0.3, 0.45, 0.6, 0.8, 1.0, 1.3, F6,_2.0, 2.5, 
3.0, 3.5, 4.0. k 3 = 0.4, \a\ = 1, ff = h=k 2 ,a k - 0, 
a 2 > 0. 
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Effect of Mach number on (a) the unsteady lift and (b) the 
unsteady moment of a Joukowski airfoil in a three-dimensional 
gust, a = 0°, thickness ratio = .12, camber = .05. fcj = 0.0, 
0.01, 0.03, 0.06, 0.1, 0.2, 0.3, 0.45, 0.6, 0.8, 1.0, 1.3, 1.6, 2.0, 2.5, 
3.0, 3.5, 4.0. h = 0.4, \a\ = 1 , ^ = - 7 , Ar x = k 2 , a ■ k = 0, 
a 2 > 0. 
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Figure 134. Effect of Mach number on (a) the unsteady lift and (b) the 

unsteady moment of a Joukowski airfoil in a three-dimensional 
gust. o = 2°, thickness ratio = .12, camber = .05. ki = 0.0, 
0.01, 0.03, 0.06, 0.1, 0.2, 0.3, 0.45, 0.6, 0.8, 1.0, 1.3, 1.6, ^2.0, 2.5 
3.0, 3.5, 4.0. k 3 = 0.4, \a\ = 1, = - j, k, = k 2 , a ■ k = 0, 

a 2 > 0. 
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Effect of Mach number on (a) the unsteady lift and (b) the 
unsteady moment of a Joukowski airfoil in a three-dimensional 
gust. a = 0°, thickness ratio = .12, camber = .10. k\ = 0.0, 
0.01, 0.03, 0.06, 0.1, 0.2, 0.3, 0.45, 0.6, 0.8, 1.0, 1.3, 1.6, 2.0, 2.5 
3.0, 3.5, 4.0. Jfc 3 = 0.4, \a\ = 1, £a. = —I, Jfc, = Jt 2 , a ■ k = 0, 
02 > 0. 



0 .25 


ct 

% 0.00 


0.25 



REAL LIFT 


0.50 


0.25 


0 .00 


0.25 



REAL MOMENT 


Effect of Mach number on (a) the unsteady lift and (b) the 
unsteady moment of a Joukowski airfoil in a three-dimensional 
gust, a = 3°, thickness ratio = .12, camber = 0. k i = 0.0, 
0.01, 0.03, 0.06, 0.1, 0.2, 0.3, 0.45, 0.6, 0.8, 1.0, 1.3, 1.6, _2.0, 2.5, 
3.0, 3.5, 4.0. k 3 = 0.4, \a\ = 1, ^ = -j, k x =k 2 ,a k = 0, 
a ■> > 0 . 
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